(N 



UNIVERSITA DEGLI STUDI DI MILANO 

Facolta di Scienze Matematiche, Fisiche e Natural! 



^ UNIVERSITA CATTOLICA DEL SACRO CUORE 



Facolta di Scienze Matematiche, Fisiche e Naturali 



O 

O 

o 

cn DOTTORATO DI RICERCA IN FISICA, ASTROFISICA E 

, — , FISICA APPLICATA 

^ Electronic structure of Ti02 thin films and 

g LaAlOs-SrTiOs heterostructures: the role of 

^ titanium 3d^ states in magnetic and transport 

g properties 

^ Settore scientifico disciplinare FIS/03 

O 
O 

o 

Coordinatore: Prof. Marco Bersanelli 

d 

(N 



Tutore: Prof. Luigi Sangaletti 



^ Tesi di Dottorato di: 

Giovanni Drera 

X Ciclo XXIV 



Anno Accademico 2010-2011 



to Luisa, Francesco and Martina 
my Best PhD achievements. 



Contents 



1 Introduction [5] 

1.1 The Ti02-5 magnetism [6] 

1.2 The conductivity of the LaAlOs-SrTiOs interface [7] 

1.3 Thesis outhne M 

1.4 Pubhcations HDl 

1.5 Acknowledgements [12] 

2 Experimental and theoretical details 1131 

2.1 Introduction [13] 

2.2 XPS and AR-XPS M 

2.2.1 Experimental set-up [15] 

2.2.2 XPS theory M 

2.2.3 Angle resolved XPS M 

2.3 XAS M 

2.3.1 Experimental set-up [23] 

2.3.2 XAS theory [25] 

2.4 Multiplet effects on core-levels spectroscopies [21] 

2.5 Resonant Photoemission [29] 

2.6 Ab-initio electronic structure calculations [36] 

2.6.1 Hartree-Fock method [37] 

2.6.2 Density Functional Theory [3H] 

2.6.3 Computing details [H] 

3 Ti02-5 Resonant photoemission [H] 

3.1 Introduction [15] 

3.2 Resonant Photoemission on Ti02 [13 

3.3 Theoretical calculations for ResPES [17] 

3.4 Experimental details and results [50] 



4 CONTENTS 



4 Magnetism in TiOa-s and NiTiOz^a films [59] 

4.1 Introduction [59] 

4.2 Growth method and characterization [63] 

4.3 Magnetic characterization [M] 

4.4 XPS measurements [7T] 

4.5 ResPES on TiOa-s and N-doped TiOa-^ 

4.6 Theoretical models [82] 

4.6.1 DFT on TiOs supercells M 

4.7 Conclusions [HZ] 

5 The LaAlOa-SrTiOa heterostructure Ml 

5.1 Introduction [HI] 

5.2 XPS data M 

5.3 XAS and XLD EM 

5.4 ResPES HUT] 

5.4.1 First data set HUH] 

5.4.2 Second data set Illll 

5.5 AR-XPS measurements 11211 

5.6 Conclusions 11321 



Bibliography 11361 



Chapter 1 
Introduction 



The transition metal oxides (TMO®) play a fundamental role in a wide 
range of technological applications. They are used as dielectrics, gate insu- 
lators and for magnetic applications. The transition metal oxides surfaces® 
have a primary importance in catalytic processes, in gas-sensing, in fuel-cells 
and in solar panels. Their conduction properties"^ range from insulators to 
semiconductors and metals, sometimes in the same compounds. This is the 
case of the so called Metal-to-Insulator transition, in which a material can 
change from an insulator to a metal state upon suitable doping or under the 
application of an external pressure. 

In each and every of these cases, much of the relevant physics is due to 
the complexity of the TM electronic structure, which still poses challenging 
problems to the physics community. In fact, while the electronic structures 
of elemental TM's have already been fairly described, the same cannot be 
stated for the TM oxides. The electron-electron interaction in the TM partly 
filled 3d shell can induce complex correlation effects, hardly described with 
the simple theories of the early solid state physics. For instance, the best 
known empirical Hamiltonian that should describe the correlation aspect of 
these systems, the Hubbard model, cannot still be solved exactly though it 
has been introduced 50 years ago. 

The complexity of the TM-oxides resides also in their bulk and surface 
crystal structure. Planar dislocations, point defects, cationic or oxygen va- 
cancies are often present in TMO crystals, but their positions in the lattice 
can be hardly controlled. The theoretical calculations on defect sites in TMO 
are also rather complex, due to the aforementioned correlation effects. Much 
of the effort of experimental and theoretical physicists (and chemists) is thus 
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devoted to the characterization and to the control of these defects. 

The present thesis is focussed on the Ti-0 bond in Ti02 based systems, 
namely Ti02-5 thin films and LaAlOs-SrTiOs heterostructures. The atten- 
tion on these system has been drawn by the occurrence of unexpected phys- 
ical properties in formally Sd" insulators: magnetism and 2D conductivity, 
respectively. Both properties are thought to be related to the 3d level pop- 
ulation of titanium cations. In the following, a brief introduction to these 
topics is provided. 

1.1 The Ti02-(5 magnetism 

A peculiar case treated in this Thesis work, is the titanium dioxide (Ti02). 
This system is a formally simple d° oxide, a transparent insulator (the gap is 
nearly 3 eV) widely applied in pigments, cosmetics and coatings because it is 
cheap, highly refractive and bio-compatible. Due to the vastity of published 
literature and to the easiness in the growth, its surface is considered as the 
prototype® of the metal-dioxide systems. The various surface terminations 
of Ti02 single-crystal are very important for catalysis applications and are 
still subject of an intensive research work. 

Titanium dioxide has also been claimed to be a room-temperature ferro- 
magnet upon transition- metal doping®. Even if at the beginning the dopant 
ions were considered as the source of the ferromagnetism (since pristine Ti02 
is diamagnetic), the researcher shifted their attention to the electronic states 
generated by the oxygen vacancies. In fact, the same kind of (weak) mag- 
netism have been recently measured in formally undoped samples; this phe- 
nomenon is known as d*^ magnetism® and has been found also in Hf02®, 
Sn02, ZnO, MgO and MgAl204®. Each of these materials, when grown with 
state-of-the-art methods, are diamagnetic. The introduction of defects that 
occurs in non-equilibrium growth techniques seems to induce the observed 
ferromagnetism. 

However, the claim of intrinsic room-temperature ferromagnetism in these 
closed-shell oxides has not yet been universally accepted by the researcher 
community®, with the most important objection being related to sample 
quality (i.e. to possible external contaminations). A fine control of growth 
and characterization is then necessary when dealing with the magnetism of 
these materials. 

The defect-related electronic states in these oxides are usually too weak 
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to be probed with standard X-ray photoemission (XPS). One of the ex- 
perimental techniques that can properly address this problem is Resonant 
Photoemission spectroscopy (ResPES), in which the valence band spectra 
are collected by tuning the photon energy close to a selected atomic shell 
absorption edge. In the case of Ti02, ResPES can enhance the titanium 
electronic states in the valence band and, if present, also in the band gap. 

In the first part of this Thesis, a comprehensive magnetic characterization 
in a set of Ti02 samples is given, together with the analysis of Ti 3d-related 
states. A series of Ti02-<5 films have been grown with the RF-sputtering tech- 
nique on various non-magnetic substrates; subsequently, a protocol of anneal- 
ing procedures has been applied in order to alter the sample stoichiometry. A 
full magnetic, crystalline and electronic structure characterization has been 
carried out at each step of the films treatments. A set of N-doped Ti02-5 
thin films have also been grown, in order to verify the effect of the N dop- 
ing in the Ti02 magnetism. The electronic states related to the presence of 
defects in these films have been further investigated through ResPES. The 
hypothesis of a clustered Vo origin of the FM is further discussed in the light 
of the experimental and theoretical results. 

1.2 The conductivity of the LaAlOs-SrTiOs 
interface 

Another interesting oxide system in which the stoichiometry of Ti ions 
play a fundamental role is the LaAlOs-SrTiOs interface (LAO-STO in short). 
LAO and STO, separately, are two band insulators, with an empty shell 
electronic structure (Sd" for STO, 4f° for LAO) and a similar perovskite 
structure; however, the interface created by growing LAO on the top of STO 
(001) has found to become metallic, hosting a quasi-2D electron gas'^^°*. This 
heterostructure becomes conductive when the STO is terminated with the 
Ti02 plane while remain insulating with the SrO plane termination; the Ti- 
related electronic states are thus expected to carry the metallic states. 

The transition from the insulating to the metallic state is observed only 
for a LAO capping larger than 4 unit cells, making the puzzle of the conduc- 
tivity even more complicated. Furthermore, in some cases a low Tc super- 
conductivity and isolated magnetic moments have also been detected, but 
not at the same time in a sample. As in the case of the Ti02 magnetism, the 
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growth conditions and the oxygen stoichiometry are crucial in the physics of 
the LAO-STO system''^. In fact the interfaces grown in a low partial pres- 
sure of oxygen yield a 3D conductivity, induced by the presence of oxygen 
vacancies. Besides, STO itself is an interesting case of defect-induced physi- 
cal properties: superconductivity ''^^ has been observed in defected SrTiOs-^, 
obtained by the annealing of pristine SrTiOs. The possibility of growing high 
quality samples has been related to the detection of the anomalous Shubnikov 
de Haas effect ^ in these interfaces, which is a clear proof of the 2D nature 
of the conductive layer. 

Different mechanisms have been invoked in order to explain the physical 
phenomena in the LAO-STO heterostructures. By recognizing that the polar 
nature of LAO could induce a diverging electric potential, an electronic re- 
construction mechanism has been proposed, based on a charge-transfer from 
the surface to the Ti atoms at the interface''^. However, there are several 
proof of cationic disorder, i.e. of a non- abrupt LAO-STO interface ^3; the 
polar discontinuity could then be energetically relaxed by the disorder and 
not by the electronic reconstruction. Finally, also the Vo's should be taken 
into account, as well as the possible charge transfer effects from surface de- 
fects''^. Each of these mechanisms don't exclude the others; at the moment 
there is not a general consensus in the researcher community. 

As in the case of Ti02, a direct experimental proof of the conductive 
electronic states can be hardly obtained with conventional photoemission 
techniques. The main difficulty is the surface sensitivity of the soft X-ray 
photoemission, since the relevant electronic states are buried beneath the 
LAO capping. Moreover, even when hard X-ray photoemission is used, the 
titanium contribution to valence band spectra is usually overwhelmed by the 
oxygen 2p states. In a similar way of d° ferromagnetism, the resonant pho- 
toemission taken at Ti edge could be the best technique to directly measure 
the Ti contribution to the conductivity''^. 

The second part of this Thesis work is devoted to the characterization of 
conductive and insulating LAO-STO interfaces, carried on with X-ray pho- 
toemission (XPS), X-ray absorption (XAS) and with ResPES techniques. 
The stoichiometry of each atomic species has been evaluated by compari- 
son with LAO and STO single crystals. A resonance enhancement of the 
conductive Ti states, associated to a small fraction of Ti'^"'" ions is reported 
and compared to theoretical calculations. In addition, a characterization of 
the intermixing and the disorder at the LAO-STO interface has been done 
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through angle-resolved XPS, which in this case has required a careful control 
over the theoretical model and over the experimental data handling. The re- 
sults are then compared to the theoretical models proposed in the literature. 



1.3 Thesis outline 

In Chapter 2 the experimental and the computational details that con- 
cern the X-ray spectroscopies of interest are presented. When possible, the 
XPS, XAS and ResPES spectra have been compared with theoretical calcu- 
lations; in this Chapter a brief introduction of the multiplet theory and the 
ab-initio density functional theory (DFT) is thus given. The Hubbard cor- 
rection to these models, the charge-transfer multiplet and the DFT+U, have 
been adopted to calculate the spectral weight and to interpret the experi- 
mental results. Finally, a section of this Chapter is dedicated to the theory 
of the AR-XPS spectroscopy, apphed to the LAO-STO interface. 

A detailed description of the ResPES spectra at Ti L2,3-edge in pure Ti02 
is given separately in Chapter 3; these data can be used as a reference 
for the experimental results given in the following chapters. The XAS and 
ResPES results have been interpreted both with multiplet and with DFT-I-U 
calculations. 

In Chapter 4, the photoemission and magnetization data for a series of 
Ti02 and nitrogen doped Ti02 films are shown. The samples have been thor- 
oughly characterized with X-ray diffraction, Raman spectroscopy, SQUID 
magnetometer and XPS at each step of annealing treatments aimed to con- 
trol the density of defects in the films. XAS and ResPES have also been 
performed on the thinnest Ti02-5 and N-doped Ti02-5 films. The role of 
nitrogen in magnetism is further explained with spin-resolved DFT calcula- 
tions. 

Finally in Chapter 5 the XPS, XAS and ResPES data taken on the 
LAO-STO samples are shown, compared to the single-crystal LaAlOs and 
SrTiOa. The experimental results of insulating and conductive samples are 
compared. The angle-resolved X-ray photoelectron spectroscopy (AR-XPS) 
and X-ray photoelectron diffraction (XPD) techniques have been used to 
check the structural properties of the interface and the capping layer. 
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Chapter 2 

Experimental and theoretical 
details 

2.1 Introduction 

A relevant part of this thesis work have been devoted to explore in details 
some core-level based X-ray spectroscopic techniques. The main advantage 
of using X-rays is the possibility to be chemically selective and thus have 
access to the electronic structure of a single ionic species. The experimental 
methods applied in this work arc X-ray photoemission (XPS), absorption 
(XAS) and photoemission taken at resonance condition (ResPES). In the 
Ti02-5 case, the magnetic characterization has been done with a SQUID 
magnetometer at Univcrsita di Pavia, as well as the Raman spectroscopy. 

X-ray photoemission measurements have been carried out mostly in the 
Surface Science and Spectroscopy Lab at Univcrsita Cattolica del Sacro 
Cuore, while X-ray absorption and resonant photoemission data have been 
collected at BACH and ALOISA beam lines at the Elettra synchrotron in 
Trieste (Italy). Electron spectroscopy spectra have been calculated with 
various electronic stiTicture calculation methods; in particular, an atomic 
multiplet approach has been adopted for transition metal core levels and an 
ab-initio density functional theory approach (DFT) for calculation of valence 
and conduction bands. The intensities of the angle-resolved photoemission 
(AR-XPS) peaks on LAO-STO have been simulated with a theoretical model 
which accounts for inelastic as well as elastic scattering. 

In this Chapter, a description of the experimental techniques and theoret- 
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ical calculations is given; some XPS and XAS results measured on manganese 
ions in germanium and cadmium telluride, already published by our research 
group, are provided as examples. 



2.2 XPS and AR-XPS 

X-ray photoemission spectroscopy is a well established photon-in electron- 
out technique that gives access to the electronic and chemical properties of 
selected atomic species on a sample surface. A schematic view of XPS is 
shown in Fig. 2.1[ b). In a XPS experiment the specimen is exposed to a 



soft X-ray source (usually hu > 100 eV) in order to induce the photoelectric 
emission from core-levels and valence band. The number of electrons vs the 
kinetic energy (E^,) spectrum is measured by an electron spectrometer. The 
survey spectra are usually composed by a series of peaks superimposed to a 



stair-like structure (see for example, the STO (001) survey in Fig. 2.2); the 
peak energy is related to a specific core level of an atoms or to secondary 
Auger electrons induced by the core-hole, while the step-like background is 
generated by the inelastic scattering of photoelectrons. XPS peaks ascribed 
to a core-level with a symmetry different than spherical s, can be split by 
spin-orbit interaction, which is inversely proportional to n quantum number; 
for example the p levels are split in pi/2 and P3/2 components, with an area 
ratio given by the electron degeneracy (1:2 for p levels, 2:3 for d levels and 
so on). 

Since the core-level energies are mostly determined by the atomic species, 
XPS allows the identification of the sample composition (chemical elements 
and their relative ratio); moreover, the specific valence state of the atom 
affects the position of the core-level peaks energies, allowing the identification 
also of the bond ionicity (or covalency) degree. This is the reason why XPS is 
also known as ESCA, Electron Spectroscopy for Chemical Analysis. While a 
complete theoretical description of the XPS technique lies outside the purpose 
of this Thesis work, some basic details about the angular dependence of the 
probing depth must be explained in order to properly introduce the AR- 
XPS method. In the next section, technical details about the experimental 
instrumentation are given. 



2.2 XPS and AR-XPS 
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Figure 2.1: Schematic view of X-ray absorption (a) and photoemission (b) spec- 
troscopies. 

2.2.1 Experimental set-up 

A standard XPS set-up requires an experimental chamber operated in 
UHV conditions (10"^° mbar) in order to avoid sample contamination and, 
to lesser extent, to avoid the photoelectrons scattering through the analyzer. 
The commercial X-Ray sources are usually based on the emission spectra of 
aluminum (Al line, hu = 1486.7 eV) and magnesium (Mg k^, hu = 1256.6 
eV), but of course XPS can be performed with the synchrotron radiation in 
dedicated beamlines, in order to have tunable photon energies and a higher 
photon flux. 

In the Universita Cattolica Labs a dual anode (Mg/Al) PSP X-ray source 
with a Scienta R-3000 AR-XPS analyzer with a 2D phosphor detector have 
been used. The best resolution was 0.7±0.1 eV, mainly due to the X-ray 
source. This analyzer can operate both in a "transmission" mode and in an 
"angular mode" that allows the simultaneous collection of the XPS spectra in 
a ±10° range. When an higher resolution was needed, a parallel XPS setup 
has been used with a monochromatized Al k^ source and a VG Mkll system, 
leading to a 0.5 eV total resolution at the expense of the total electron count 
rate. 
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Figure 2.2: XPS survey spectrum on the STO (001) surface, taken with a 
ka source. 



2.2.2 XPS theory 

From the conservation of total energy, the kinetic energy of an electron 
photoemitted from a core-level can be expressed as follows: 



E 



Kin 



hv - Esin - ^ 



(2.1) 



where E^j^ is the electron binding energy and $ is the work function, which 
is the extra energy needed to transport an electron from the sample to the an- 
alyzer. The electrons from Auger process are related only to internal atomic 
relaxation and thus are photon-energy independent (i.e. fixed kinetic en- 
ergy), while the kinetic energy of core-level photoelectrons changes with the 
photon (fixed binding energy). A tunable (multiple anode, or synchrotron 
based) X-ray source can thus be used to better identify Auger structures, 
when superimposed to other core-levels. 

The probing depth of this technique is due not to X-ray penetration (in 
the ixm range) but to the inelastic mean free path (IMFP) of electrons in 
solids which, at these kinetic energies, is in the order of 1-2 nm: XPS is 



thus a surface-sensitive technique. In Fig. 2^ in section 5.4 are shown some 
calculated IMFP for the LaAlOa case. In some synchrotron facilities, the 
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photoemission with hard X-ray photons (HAX-PES, hu up to 10 keV) is also 
possible, which in turn is a bulk- sensitive technique; in HAXPES however 
one has to deal with a reduced X-ray flux and, most important, with smaller 
photoelectron cross-sections. 

From a quantum-mechanics point of view, the photoemission process can 
be written as follows with the Fermi Golden rule: 



Where the squared matrix element gives the transition rate and the Dirac 
delta accounts for energy conservation; in photoemission, the flnal state is 
the ground state plus a core hole and a free (photoemitted) electron with 
Efc energy. The matrix operator T for the electron-photon interaction can 
be taken as the simple dipole operator T = e ■ r , since the contribution 
of the quadrupole operator for the soft X-ray regime is neglig ibleCS. This 
equation can be solved with the Green functions method, in order to account 
for scattering process and flnal state effects due to the relaxation of the 
electronic levels in the proximity of the core-hole; this approach is called the 
"one-step" photoemission theory. 

Another way to describe XPS is the so-called "three-step" model, in which 
the process is described in three different steps: 

• flrst step is the photoemission from the atom, described with photoe- 
mission cross section; 

• the second step is the drift of the electron in the solid, described by the 
electronic IMFP; 

• the third step is the escape of the photoelectron from the solid, de- 
scribed by the work-function. 

The three-step model is often assumed in the framework of the sudden 
approximation, in which one suppose that the core-hole flnal state doesn't in- 
fluence the XPS spectra, while the one-step model gives e a good description 
of excitonic effects and plasmon resonances. 

In the three-step model, the intensity (i.e. the area) of a photoelectron 
peaks depends on many parameter; disregarding the diffractive effects (very 




(2.2) 
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important in single-crystal, as will be shown in Chapter 5) and X-ray attenu- 
ation (important only for grazing photon incidence angles), the contribution 
to the photoemission intensity of an infinitesimal thick layer at depth z from 
the surface can be expressed as follows: 

dl = (f)Do{EK)(T{hu)pAoP{X, z, 9)dz (2.3) 

where 

P{X,z,e) ^e~^^ (2.4) 

and: 

• is the photon flux; 

• Dq is the analyzer transmission function that depends on the kinetic 
energy (old electron analyzer could have different sensitivity at different 
Efc); 

• (J is the photoelectron cross section that depends on the atomic species, 
the electronic level (i.e. ls,2s,2p...) and the photon energy; 

• pAo is the number of atoms per unit volume per sampling area, that can 
depend on the X-ray focalization on the sample and/or on the analyzer 
focus; 

• P{X,z, 9) is the probabihty of an electron at depth z to escape from the 
specimen and to reach the analyzer placed at an angle 9 relative to the 
sample normal; A is the IMFP and depends on the sample composition 
and on the E^.. 

When dealing with the elements ratios only, the photon flux and sampling 
area can be discarded since they are contributing in the same way to each 
photoelectron peak. When measuring XPS peaks with close kinetic energy 
(or with a well calibrated analyzer) also the Dq term can be neglected. The 
other parameters have to be taken into account to carefully quantify the ele- 
ments concentration; the typical accuracy of this method in a homogeneous 
sample is around 1-5%. The typical XPS low boundary of sensitivity is about 
0.1-1%. 
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2.2.3 Angle resolved XPS 



The parameter P in Eq. |2.4| is also dependent on the relative orientation 
of specimen with respect to the analyzer; in fact, at grazing emission, photo- 
electrons of inner layer have to travel for a longer path through the sample 
in order to reach the vacuum, resulting in a reduced photoelectron intensity 
from the bulk. The surface sensitivity of XPS is thus improved at grazing 
emission angles; this can be useful for instance for enhancing the contribution 
of absorbed molecule respect to the host substrate. 

This effect can be also used in a reverse way: for example, in the case of a 
thin film coverage of a bulk material, one can use AR-XPS data to estimate 
the overlayer thickness. This is only possible when the composition of each 
layer is known, in order to fix all the parameter described in Eq. |2.3t in some 
cases just a single measure (i.e. taken at only one emission angle) could be 
sufficient to estimate thickness, but of course a set of data taken at different 
angles can improve the results reliability. 



Eq. 2.3 has to be integrated on the z axis to obtain the total XPS in- 



tensity; for instance, in the case of a multilayer sample on a substrate (see 



Fig. 2.2.3) with thickness di and the contribution on the photoelectron 



intensity of each layer can be written as follows: 



1 P . cosid) , , , 

ll^aaPa e — wd2; = aapa— ^(1-e '^i-w) (2.5) 







2 f'^'-^'^^ ^ , Cos(9) A, 

I^(X(Tbpb / e -cos(o) = abpb — -^e diCo.m (^i _ g daCo^w) (2.6) 



Ac 



di+d2 -^c 

where a, b and c are index for atomic levels; the photon flux has been omitted, 
since only the relative ratios are usually evaluated. These equations can be 
used as a model to flt AR-XPS data in order to obtain the di and ^2- With 
simple modiflcation of the element proflle along the z axis one can deal also 
with other cases, such as 3D islands on a surface. The photoelectron cross- 
sections in the soft X-ray regimes for each element of the periodic table have 
been calculated by Yeh and Lindau'^; the IMFP data can be evaluated 
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Figure 2.3: Schematic 
configuration of an AR-XPS 
model for a multilayer sample. 
In this Figure the parameters 
2 3|2.6|2 .7 



for Eqs. 



are defined 




through the Tanuma-Powell-Penn® (TPP-2M) formula, that requires the 
knowledge of material density, bandgap and the number of valence electrons, 
as well as the electron kinetic energy. 

Unfortunately, this model suffers from two drawbacks. First, different 
topological arrangements can lead to exactly the same peak ratios, as shown 



in Fig. 2.2.3 a comprehensive approach should consider also the slope of the 
stair-like background, which increases with the number of inelastic scattering 
events®. 

Besides, even if the structure composition (overlayer, multilayer, 3d island 
etc.) were known, the inelastic mean free path alone would not be sufficient 
to correctly estimate the thickness of the multilayer structure. A soultion 
to this problem is given by the replacement of IMFP with the Effective 
Attenuation Length ''^3 (EAL), which is a more complete estimation of the 
effective opacity of a solid for a given electron energy. The EAL is calculated 
by introducing the contribution of the elastic scattering length, which can 
change significantly both the exponential and the angular attenuation of the 
XPS signal. As a matter of fact, the real attenuation length is dependent to 
the specific geometry of the system and thus can be only evaluated through 
computational demanding Monte-Carlo calculations. However, it is possible 
to obtain EAL values for specific cases (like thin overlayer or thin marker- 
layer at a specific depth) through simple formulas. EAL values adopted in 
this Thesis have been calculated with the EAL program of NIST''^^. The 
expression for the local EAL are: 



EAL = ~[cos{a) 



dln(j){z, a) 
dz 



(2.8) 
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Figure 2.4: Different sur- 
face structure of copper in 
gold which give exactly the 
same XPS Cu 2p peak inten- 
sity. The different configura- 
tions can be inferred by inelas- 
tic background at higher BE 
than the photoemission peak. 
From Ref . CD 



for measurement of marker-layer depth and 



d f°° 

EAL = —[cos{a)—ln( / (f){z,a)dz)]~'^ 
at 



(2.9) 



for measurement of an overlayer-film thickness t, where (j) is the emission 
depth distribution function that depends on depth z and on the analyzer 



angle a. The EAL from Eq. 2.9 can be used instead of A in formulas such 
as Eq. 2.7 before the integration, which thus can't be done analytically (the 
EAL depends on z). An average value of the EAL, designed to be used for 
convenience in the XPS intensity formulas after the integration, can also be 
obtained through the EAL code; however, this "practical EAL" is suitable for 
a thick (more than 3 nm) capping layer and thus has not been adopted in 
this Thesis work. 

As an example, in Fig. 2.5 the graph of EAL vs depth for electrons 
photoemitted from La 4d and Al 2s core levels in LaAlOs is shown. EAL is 
usually lower than the corresponding IMFP and rapidly decreases in a short 
distance from the surface. As it is shown in Chapter 5, the introduction of 
the EAL is thus fundamental to obtain the correct peak intensity or, on the 
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Figure 2.5: Effective attenuation lengtli jgp^]^ (i.e. tfie distance from 
the surface) calculated for La 4d and Al 2s core levels photoelectrons on LaAlOs; 
the dotted lines indicate the corresponding inelastic mean free path, calculated 
with TPP-2MCSJ formula. The calculations have been done for an Al ko. X-ray 
source. 



contrary, the thickness of a very thin overlayer. 



2.3 XAS 

In a XAS experiment the X-ray total absorption cross-section of the spec- 
imen is measured. Therefore, XAS requires the possibility to scan the photon 
energy and thus it is usually performed on synchrotron facilities, even if sim- 
ilar spectra can be obtained through the electron energy loss spectroscopy 
(EELS) which formally share the same theoretical formalism. While XAS is 
the most generic name, many different notation are given according to the 
required experimental information or to the specific application fields; here 
is a list of the different XAS definitions: 

• NEXAFS (Near-Edge X-ray Absorption Fine Structure) or XANES 
(X-ray Absorption Near-Edge Structure): in both cases the absorp- 
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tion cross section is measured "near" the absorption edges of a specific 
elements. Altfiough from an experimental point of view are both sim- 
ilar, in general NEXAFS is used in surface or molecular studies and 
XANES in crystal or bulk studies. The NEXAFS experimental infor- 
mation covers the empty states, the adsorbate geometry (NEXAFS) 
and the symmetry of a specific atomic species (XANES). 

• EXAFS (Extended X-ray Absorption Fine Structure): in this case, 
the absorption cross section is collected over a wider energy range. 
In ordered structures, these spectra show oscillations due to multiple 

scattering of excited electrons: by analyzing the Fourier transformation 
of the EXAFS spectra it is possible to evaluate the nearest-neighbors 
bonding lengths. 

• SEXAFS (Surface Extended X-ray Absorption Fine Structure): same 
as EXAFS, but tuned to give more information on the surface bond 
lengths. 

• XMCD and XLD (X-ray Magnetic Circular Dichroism and X-ray Lin- 
ear Dichroism): XAS spectra taken with different X-ray polarization. 
XMCD requires a magnetic field (both of specimen or external) and 
gives direct information on the local magnetism; XLD can give fur- 
ther information on the symmetry distortion of the ionic environment. 
These techniques can be done also in an imaging-hke fashion, both 
scanning the sample position or using electronic lenses, giving the pos- 
sibility to map the specimen different magnetic or structural domains. 

In this Thesis XANES, XMCD and XLD measurement are reported. Fur- 
ther details about experimental set-up and theory will be given in the next 
paragraphs. 

2.3.1 Experimental set-up 

Most of XAS data shown on this thesis have been taken on BACH beam- 
line at ELETTRA synchrotron in Trieste (Italy). The reference spectra taken 
on Mn:CdTe and Mn:Ge have been measured at the ALOISA beamlines, 
again at ELETTRA. A synchrotron X-ray source combines an high brilliance, 
a small on-target X-ray focus as well as the photon-energy and polarization 
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Figure 2.6: Top view of BACH beamline layout at ELETTRA synchrotron. From 
left to right are depicted the undulators (Ui and U2) on the storage rings, then 
the mirrors, the slits, the monochromator and the user end-stations. 



tunability. In Fig. 2.6 a schematic view of the BACH (Beamline for Ad- 
vanced diCHroism) beamline set-upis given. The photon energy is tuned by 
changing the distance (usually referred as "gap" ) of undulators magnets and 
by means of a two mirror monochromator. The resolving power (AE/E) at 
BACH beamline is 20000-6000, 20000-6000, and 15000-5000 in the energy 
ranges 40-200 eV, 200-500 eV, and 500-1600 eV, respectively. As in XPS, 
XAS requires UHV conditions and, because of synchrotron high photon-flux, 
can suffer of charging effects. The X-ray polarization can be tuned with 
dedicated insertion devices (helical undulators) or by exploiting the natural 
polarization characteristic of synchrotron light. 

The absorption cross section could be in principle detected by measuring 
the photon flux before and after the sample; in practice, this method can 
be applied only to very thin (in the range on /im) samples. Usually XAS 
is performed by measuring secondary de-excitation process caused by the 
absorption of X-rays. Here is a summary of detection technique: 

• Fluorescence: XAS can be measured by detecting the rate of fluores- 
cence given by the recombination of the electrons with the core-hole 
created by the excitation. This detection technique requires the pres- 
ence of a silicon-based photon detector in the measurement chamber 
and its sensitivity is related to the X-rays penetration depth in the sam- 
ple (in the /im range). Fluorescence detection is less effective in light 
materials, since the Auger decay is the most probable de-excitation 
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process, and can be quenched in very dense material, because of the re- 
absorption of the emitted photon. It doesn't suffer of charging effects, 
even in insulating samples. 

• Total yeld: XAS can be detected by measuring the electrical current 
( "drain current" ) generated by the X-ray absorption. A picoammeter is 
needed, since this current ranges typically on the 10^^*^- 10^^ A scale. 
The drain current is generated by a cascade of Auger process that 
are also related to the electrons inelastic scattering. Since only the 
electrons that reach the surface contribute to this current, the probing 
depth is lower than in fluorescence detection. The typical total yield 
MPD (maximum probing depth) is in the 4-10 nm range, the Ti L-edge 
being at the lower MPD side''^. 

• Partial yeld: the photoemission intensity is usually proportional to the 
absorption cross section, thus XAS can be measured by integrating the 
photoelectron emission in a defined energy range, through a channel- 
tron detector or an electron analyzer. This detection technique is the 
most surface-sensitive and, like XPS, can be affected by charging effects 
in insulating samples. 

Spectra should be normalized with the incoming photon flux, which is 
usually measured on the monochromator last mirror through the drain cur- 
rent method. 

2.3.2 XAS theory 

A schematic view of XAS is shown in Fig. |2.1[ a). In short, when the 
photon energy is higher than a core-level binding energy, an electron from 
that core level could be excited into an empty states below the Fermi edge. 
The transition rate for this process can be described with Fermi Golden rule, 
similarly to Eq. |2.2 

WpU = ^mT\f)\^ p{Ef-E,-hu) (2.10) 

where p is the empty level density of states (DOS). As in XPS, operator T can 
be taken as the usual dipole operator e-r. In the case of most core- level edges, 
except in transition metals and rare earths, with photon energy between 
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Figure 2.7: Rutile Ti02 K-edge XAS spectrum, as compared to ab-initio DFT 
calculations. In blue, the calculated p-projected density of states is shown. Taken 
from Ref.CS. 



100-2000 eV the dipole operator is slowly varying with hu; in this situation, 
the XAS spectra becomes a direct measurement of the empty DOS, plus a 
contribution from electrons multiple scattering . A proof of this can be 
seen in Fig. 2.7, where an oxygen K-edge XAS spectrum on Ti02 is compared 
to the DOS results of a DFT ab-initio calculation. The dipole operator is still 
important, since it gives the selection rules and the polarization dependence. 

In the case of open-shell system, such as transition-metals (TM) or rare- 
earths (RE), the electrons excited in the empty states can strongly interact 
together and with the core-hole. In such system, the matrix element of 



Eq. 2.10 becomes the leading term and multiplet features appears in the 



spectrum. Since many experimental data of this Thesis have been taken at 
the titanium 2p threshold, the multiplet effects are explained in more details 
in the next section. 



2.4 Multiplet effects on core- levels spectro- 
scopies 



There are certain cases in which the XPS and XAS matrix elements of 
Eqs. 2.2 and 2.10 have a complex photon energy dependence; in fact, in the 
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case of TM and RE multiple final states are allowed, because of the electron- 
electron interactions in the empty d shell (for TM) or f shell (for RE). In the 
TM cases the 2p, 3p and even the 3s level XPS and XAS spectra are thus a 
superposition of many electronic transition between the core-hole level and 



the 3d shell (see Fig. 2.8) 



Simple ground-state DFT calculations cannot account for these effects be- 
cause they are deeply related to the excited states; on the contrary, Hartree- 
Fock (HF) calculations can be done also for excited states, even if the elec- 
tronic correlation is only partly described. The latter approach (HF), mixed 
with a proper electronic angular momenta algebra''^, has proved to be an 
effective tool in predicting core-level spectroscopic data, under the name of 
atomic multiplet theory. This model can be expanded with subsequent steps 
in order to better match the experimental results on crystals. 

The starting point is an HF calculation for a single ion, carried out with 
the COWAN ^ code, which gives the energy levels and the transition ma- 
trix from the ground to the final states. The results is an atomic multiplet 
spectrum that describes the XAS (or XPS) process in a spherical symmetry 
(SO3). The eigenvalues and the matrix elements are then modified according 
to crystal field (CF) theory and an finally to an Hubbard model to mimic 
the ionic environment. Atomic multiplet plus CF calculations have been 
carried out in this Thesis with the MISSING package and the full calcu- 
lations, which also include the ligand-metal charge-transfer (CT), have been 
performed with the CTM4XAS ^ code. The complete process is shown for 



a 3d° atom in Fig. 2.8 To account for both inter-atomic charge-transfer 
and intra-atomic multiplet splitting effects, a configuration interaction (CI) 
description of the wavefunctions is usually adopted. Several configurations, 
denoted as rf", d"~^^ L, c?""''^/^ and so on {L denotes a ligand hole) are 
used to describe the open shell of the 3d transition metal ion during the 
photoemission process. Accordingly, the initial state wavefunction is written 
as: 

= ai|3rf") + a2|3rf"+^L^) + a3\3d''+^ l/) (2.11) 

where L denotes a configuration with a p-hole in the anion states. In the 
present case, the p-hole represents the CT from the O 2p levels to TM 3d 
levels. The CT energy is defined as A = E{d"'~^^L) — E{d'^), whereas the 
Coulomb d-d interaction is represented hy U = E{d"-^^) + E{d"-^^) —2E{d"') , 
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where E{d'^L'^) is the center of mass of the d'^L'^ multiplet. The Coulomb 
and CT energies are usually regarded as model parameters. Interaction be- 
tween d'^^^L'^ and (^"+'"+12, "i+i configurations is accounted for by the Tp^ 
off-diagonal term in the Hamiltonian matrix. 

In the case of core level photoemission the final state interaction between 
the core hole and the 3d electrons in the outer shell is explicitly accounted 
for by an energy parameter Q. Therefore the CT term is corrected as A — Q 
upon creation of the core hole. The final state wavefunction is given by: 

= /5i|2p3rf") + /32|2p3rf"+^Li) + P^\2p?>d''+%^) + e" (2.12) 



where c represents the core hole. The spectral weight in a photoemission 
experiment is calculated, in the sudden approximation, by projecting the 
final state configurations on the the ground state, i.e. 

/xps(5^) oc5^|(^G5|^^,/s)|'5(5E-£,) (2.13) 



where 



|(*Gs|^.,/s)|' = + a2/32,i + • • f (2.14) 



and the sum is run over all final state configurations \^ijs) with energy 
Si. Each of the d" configurations can be further described with the crystal 
field theory, with the multiplet calculations or with both approaches. The 
CF is introduced by reducing the symmetry from the spherical SO3 group 
to the desired symmetry; in this thesis, an octahedral {Oh) and tetragonal 
(D4/1) space group are used. The adequate level of the theory approximations 
is fixed by the experimental technique. For example, in XAS a multiplet-CF 
approach is usually sufficient to describe most of the experimental features, 
while in XPS the CT effects cannot be neglected. 



An example of these calculations is given in Fig. 2.8 for Ti^"^ (3d°) 2p 
XAS and XPS in an octahedral symmetry {Oh}- Without multiplet, only 
two peaks separated by the spin-orbit interaction are predicted. The interac- 
tions in the 3d levels (and with the 2p core-hole) induce a complex multiplet 
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structure which is further spht by the introduction of the CF. Eventually, 
the addition of charge-transfer configurations adds small (strong) satellites 
to XPS (XAS) spectrum. As can be noted in Fig. 2.8, the shape of the Ti 
L2,3 edge XAS spectra is mostly determined by CF effect, which are less im- 
portant in XPS; on the contrary, CT satellites are best observed in XPS. The 
small CT satellites in XAS are due to the conservation of the total charge in 
the absorption process. 

The relative intensities of the XPS satellites with respect to main peaks 
are related to the type of chemical bond. In more ionic compounds (such 
as MnO) weak satellites are observed, therefore the spectra can be often 
described with multiplet calculations only; in samples with covalent bonds, 
strong XPS satellites are present and thus in order to evaluate the exper- 
imental spectra the charge-transfer model is required. As an example, in 



Fig. 2.9 multiplet CT calculations (without CF) and experimental data''^ 
for manganese ions in different host matrices are shown. The shift from an 
atomic-like spectrum (small satellites, in MnO) to a covalent bond (large 
satellites, Mn:CdTe) is clearly visible. As shown in Fig. 2.9, the charge- 
transfer parameter (A) determines the strength of the satellites peaks. In- 
sulating samples where A < U are referred as charge-transfer insulators''^; 
in that case, which includes Ti02, A is close to the band-gap (see the A 



parameter in Fig. 2.8). In the case of a metallic material the XPS peaks 



assume a typical asymmetric shape (Doniach-Sunjic''23), which cannot be 
described with atomic calculations. 

The ligand-field charge-transfer multiplet calculations provide an effec- 
tive tool for calculating XAS and XPS spectra in TM; however, this theory 
requires numerous empirical parameters: a reduction factor for the Slater in- 
tegrals (obtained by HF), the crystal field parameters related to the crystal 
point group (i.e. lODq, ds,dt and so on) and the U, A, Q an T parameters set. 
Even if most of these parameters can be obtained from various experimental 
techniques, a complete ab-initio approach would be highly desirable. 



2.5 Resonant Photoemission 

In a ResPES experiment, the valence band photoemission spectra are 
collected by changing the photon energy across an absorption edge. In a 
single-particle approach (i.e. multiplet calculations and following approxima- 
tions) the direct valence band photoemission channel (c'^w" — )■ c^w""^ + e~) 
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Figure 2.8: Step in charge-transfer multiplet calculation of 2p XPS (left) and 
XAS (right) spectra for a Ti'^''' ion. The upper graphs are the atomic multiplet 
calculation, the middle ones are the results after the introduction of an octahedral 
crystal field and the lower ones are results with the charge-transfer (CT). Black 
arrows indicate the CT satellites. CT parameters are 10Dq=1.7 eV, A=3.0 eV, 
Q=5.0 eV, U=9 eV and the mixing parameters are 3.0 and 1.0 eV for T2g and 
symmetries, respectively. 
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Figure 2.9: Comparison of Mn 2p XPS spectra from MnrCdTe (100), 
Zni_2;Mn3;Te, MnO and a reference metallic Mn thick film deposited on a sili- 
con wafer. Calculated Mn 2p XPS spectrum for Mn^+ ion with charge-transfer 
derived configurations (A=2.1 eV, T=2.2 eV, Q=6.0 eV, and U=5.1 eV). Taken 
from Ref.^23. 



can interfere with the autoionization channel caused by the presence of a 
core-hole (c'^-i|;"+i _ > c™3(i"~^ + e~). As a results of this, the part of the 
valence band related to the specific atomic species can be greatly enhanced 
(or suppressed). This effect is known to occur for most of the elements, both 
organic and inorganic. 

Following the description of Briihwiler et al.®, a pictorial representa- 
tion of the process is given in Fig. 2.10 In this figure the most important 
excitation-deexcitation processes related to the ResPES technique are shown. 

On the left part the photoemission-related channels are shown: Fig. 
2.10K a) shows the normal valence band photoemission (VPES) and Fig. 



TWjh) the core-level photoemission. In the latter case, the system can 



fill the core-hole through a normal Auger decay (Fig. 2.10[ d)) or through 
fluorescence, which is not usually measured during a conventional ResPES 
experiment. In Fig. 2.10K d) the resonant transition from the core level to 
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Figure 2.10: Schematic of the possible excitation and de-excitation channels in 
a ResPES experiment. 



the empty state is depicted (i.e., the X-ray absorption). In this case the sys- 
tem can relax by filling the core-hole with a valence band electron, opening 
two different autoionization channels: Fig. 2.10[ e) shows the emission of the 
excited electron from the empty states, labeled as partecipant decay and Fig. 
2.10[ f) shows the emission of another electron from the valence band, labeled 



spectator decay. The final state of the participant decay is equivalent to the 
valence band photoemlssion one, leading to the (constructive or destructive) 
interference effect that is called Resonant PES (RPES or ResPES). 

The spectral weight related to the spectator decay and to the normal 
Auger channels are usually similar in shape, with a shift in energy due to the 
lower screening effects induced by the extra electron in the conduction band 
(see Fig. [211] ). 

The ResPES technique can give information also on the excited electron 
dynamics, because the intensity of the interference effects between VPES 
and parteclpator channels is related to the excited electron lifetime. In fact, 
in an highly delocalized (metallic) band the extra electron can be quickly 
removed, quenching the ResPES channel to the normal Auger. ResPES 
is thus a useful experimental tool to evaluate the charge-transfer process, 
especially between a substrate and deposited molecules or atoms. When 
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Figure 2.11: The Auger (dots) and resonant valence band photoemission spectra 
(thick line) on fullerene CeO. The energy difference between Auger and spectator 
decay, called spectator shift is clearly visible (adapted from Ref. ''^). 



applied for this purpose, the ResPES technique is usually referred to as Core- 
Hole-Clock (CHC) spectroscopy. The interference intensity can be calculated 
in the Fermi golden rule approach as follows'^: 



f 



{f\Vr\9) + Yl 



{f\VA\m){m\Vr\g) 



En 



E — ?T 



/2 



KEf - E,) 



(2.15) 



where the left part of the matrix element is the transition rate for nor- 
mal photoemission from ground (g) to final state (f), the right part is the 
summation over the possible intermediate state (m) and is the lifetime 
of excited state. Vr and are the radiative (dipole) and the Coulomb 
(Auger) operators. Under simple assumptions (only one core-hole excitation 
and a continuum of empty states) this formula can be simplified in the Fano 
formalism , obtaining: 



Y.\^f\Vr\g)\' 



{q + e) 



Em/H 



l + £2 {Ef + Eg-E^Y + Vr 



(2.16) 
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Figure 2.12: Effect of Fano q-parameter on ResPES intensity. 



where the ResPES intensity is the normal VPES one multiplied by the factor 
{q + e)"^ / {1 + e"^) with e = {E — Em) /T . The parameter q is called Fano factor 
and describes an antiresonant effect for g = and an high resonant effect 
for high g > 2, as can be seen in Fig. 2.12 In the case of transition metals 
(TM) and rare earths (RE) the q factor can be rather high (g > 3), leading 
to a pronounced resonance effect called "giant resonance" . 

An example of ResPES, in the case of Mn over a CdTe (110) surface''^, 



is given in Fig. 2.5 A clear resonance effect is seen in Fig. 2.5 the on 



resonance spectra is different in shape and overall intensity from the off- 
resonance spectra. The difference between the two (spectrum (e) in Fig. 2.5) 
is related to the partial DOS of manganese 3d electrons in the valence band 
that, in this case, can be calculated with a multiplet approach. At least, two 
final state configurations have to be considered in order to reproduce the ex- 
perimental data: the expected 3d^, that results from a photoemission process 
3(i^ — )• 3(i^+e~ resonating with the autoionization process 2p^?)d^ — )■ 3(i^ + e~, 
and a charge-transfer 3d^L configuration due to the Mn-Te hybridization. In 
order to evaluate the non-resonant valence band contribution a more detailed 
method should be used, such as ab-initio DFT. 

Another interesting case where "valence band" (i.e. 3d electrons) multi- 
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Figure 2.13: Resonant photoemis- 
sion data of MniCdTe taken at Mn L2,3 
edge, (a) Two dimensional plot of Re- 
sPES data for photon energies ranging 
from 636 to 655 eV; (b) XAS spec- 
tra at Mn L2,3; (c) and (d) are the 
on-resonance and off resonance valence 
band spectra, and their difference (e) 
is fitted with multiplet calculation (f) 
composed by d^ and d^L configura- 
tions. Adapted from Ref. ® 
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Figure 2.14: Valence band spectra on MnsGes taken at Mn L3 edge, as well 
as with hz/=48 and 51 eV (i.e., across the Mn 3p-3d resonance), and calculated 
photoemission spectrum for the Mn"*" 3d^ ion in the atomic-multiplet approach. 
Adapted from Ref.®. 



plet calculations can reproduce the average energy position of the resonating 
features in the VB is the MnsGes alloy. In Fig. 2.14 are reported the res- 
onant spectra''^ taken at Mn L2,3-edge, with theoretical calculation; the 
metallic character of this material can be simulated by considering, instead 
of a Mn^"^ (Sfi^) ion, a Mn^"*" {3d^) ion without the spin-orbit interaction 
on the d-level. The same approximation has proved to be effective also in 
calculating core-level photoemission and absorption. 

Since many experimental results of this Thesis work have been achieved 
through the ResPES technique at Ti L-edge, the Chapter 2 is dedicated en- 
tirely to a detailed description of the ResPES process in the Ti02 d° system. 



2.6 Ab-initio electronic structure calculations 

As already pointed out, "simple" multiplet calculations can deal very well 
with transition matrix to excited states in open-shell system, but fail in re- 
producing completely the valence band spectra and the DOS in general. For 
this reason, ab-initio''^ calculations have to be resorted, with a complete 
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description of the wavefunctions that participate to the chemical bond. The 
problem addressed in ab-initio calculations can be described with the Hamil- 
tonian for N interacting electrons (in r position) in a lattice of M ions (in R 
position) : 

N N,M N 2 

2m ^ * ^\ri-Ri\ ^ n - r,- ^ ^ 

i i,j ' •' ' i<j ' ■' ' 



= (2.18) 

This is already an approximation (Born-Oppenheimer), since it lacks of the 
kinetic term for ions and the ion-ion repulsion (the latter can be taken into 
account with empirical formulas). However, even with these approximations, 
the solution of this Hamiltonian is a formidable task that requires intensive 
numerical calculations. Many advances in the field of quantum chemistry, 
such as DFT, the LDA approximation, pseudopotentials and so on, allow 
now to perform complex (and realistic) ab-initio calculation even on a single 
desktop/workstation machine. In this Thesis, various ab-initio calculations 
results on simple and doped crystals are shown, mostly based on the Density 
Functional Theory (DFT). A brief description of this method is given in this 
section, starting from the Hartree-Fock equations. 

2.6.1 Hartree-Fock method 



The first successful attempt to solve the Eqs. |2.18| and |2.17| was the 
Hartree-Fock (HF) method®, in which the total wave function \1' is ex- 
pressed as a Slater determinant of single particle functions: 

^ = Det\iJi{rj)\ (2.19) 

This formulation gives a correct antisymmetric description of the total wave 
function. It is possible then to define the total energy of the system and to 
minimize it through Lagrange multipliers in order to obtain the correct single- 
particle wavefunctions ip{r) and the corresponding energies. By knowing the 
total energy of the system, it is possible to find the lattice (or molecular) 
equilibrium position, the binding energies and even the phonon frequencies; 
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by knowing the single-particle wavef unctions, it is possible to compute also 
the band structure, the density of states and evaluate the stoichiometry of 
each atom. 

This method, at least in the initial formulations, suffers from some draw- 
backs: the exchange interaction energy has to be computed with a two par- 
ticle integral (since the exchange interaction is non-local), which has to be 
often approximated to a single particle one to speed up calculations; more- 
over, the HF wavefunction cannot properly describe the electronic correlation 
because of the oversimplification of a single-particle approach. In order to 
include the correlation to HF, Configuration Interaction (CI) methods have 
been introduced®, in which the total wave function is described as a sum 
of many Slater determinant. The expansion of the ground and excited state 
in charge-transfer multiplet calculations is reminiscent of this approach. The 
fuU-CI is probably the most accurate electronic structure method but it is also 
very computational demanding, and can be applied only to small molecules 
(around 10 electrons). 

2.6.2 Density Functional Theory 

Another answer to the problem was given by the Density Functional The- 
ory (see Ref. ''^ and Refs. therein), which is probably the most popular 
electronic-structure ab-initio method at the moment. The theoretical foun- 
dation of DFT are given by the Hoenberg-Kohn (H-K) theorem, which states 
that in a system of interacting electrons in an external potential (in this case, 
the potential generated by the ion lattice) there's a bijective correspondence 
between the electronic density and the external potential. An equivalent 
statement is that the ground-state energy of such systems, apart for an un- 
interesting additional constant, is a functional of the total electron density n 
solely: 



where T is the energy associated to kinetic term, U the contribution of the 
electron-electron interaction (both Coulomb and exchange) and Yext is the 
external potential. 




(2.20) 
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By means of the variational theory it should be possible to find the 
ground-state electronic density of the system by minimizing the total en- 
ergy. The power of the H-K theorem resides in the fact that the total energy 
can be computed out of a 3 variable function (the total density) instead of 
the complex 3N \l/(ri, r2, t„) wave function. However, even if V[n] is a 
simple term, it is rather difficult to express T[n] and U[n] as analytical density 



functionals, therefore the direct minimization of Eq. 2.20 cannot be easily 
carried out without further approximations. 

This equation can be solved instead in the so called Kohn-Sham (KS) 
scheme, in which the total density is expressed through a fictitious single- 
particle system in the form of Slater determinant of tpi, functions. The 
T[n] is thus broken in two parts: 



h ^ r 

T[n] = T,[n] + T,[n] = -^Yl I d'rV'^^{r) + T,[n] (2.21) 

i 

where Th is the "Hartree" part of kinetic term, while Tc is the extra correla- 
tion energy which is unknown but, due to H-K theorem, is still a functional 
of the density; the same can be done with the U[n] term, which is splitted in 
a Coulomb part (U^) and an exchange part (Ug). The total energy is then 
rearranged in this form: 

E[n] = T[n] + U[n] + V[n] = T^[n] + Uh[n] + J d''rn{r)v,,t + E,,[n] (2.22) 



where E3;c[n]=Tc[n]-|-Ue[n]. The Eq. 2.22 is formally correct, except for the 



fact that exchange-correlation energy (E^ic) is not known. The real success 
of DFT came with a suitable (and effective) approximation the exchange- 
correlation potential Vxc (i-e. 6Exc/Sn): the local density approximation 
(LDA). In this approximation, the Va;c[n(r)], functional of the position and 
of the density at that position, is equal to the exchange-correlation potential 
of an homogeneous electron gas (HEG) with the same electron density. 

Data of V^^*^ as function of the density have already been calculated 
with numerical Monte-Carlo methods by Ceperley and Alder '^9, and thus 
can be easily interpolated with a parametric function and implemented in 
the calculation software. Eq. |2.22 thus can be solved through self-consistent 
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methods like HF but in a faster way because in DFT the exchange potential 
is fully local. Moreover, the results contain (in an approximate way) also 
the correlation energy, introduced by the LDA formalism. However, when 
dealing with TM or RE, the correlation effects became very strong and further 
refinements have to be included in the theory. 

Even if LDA should be a rough approximation, since electron density is 
not a slowly varying quantity in molecules or crystals, the results in terms 
of bonding length, binding energies and other properties are usually rather 
good. DFT can be also applied to magnetic systems, by considering two 
separate spin-up and spin-down electronic densities (n^ and n-^) within the 
Local Spin Density Approximation (LSDA). Ground state DFT calculation 
are now considered as the starting point for more accurate theoretical models. 

The improvements of DFT have gone in two directions. On the first side, 
the DFT has been extended to excited states (since KS-DFT is a ground- 
state theory) and to highly correlated materials. In fact, with a correct 
description of excitations it is possible to properly evaluate the band-gap, 
the optical spectra and other relevant physical properties. A well known 
DFT formulation for excited states is the Time Dipendent DFT (TD-DFT). 
In highly correlated materials, a better agreement with the experiments can 



be achieved by adding extra-energy terms to the Eq. |2.22[ such as in the 
GW approximation (a self energy correction) or in the DFT+U (an extra 
Hubbard term in selected energy levels). 

On the other side, many attempts have been done to improve the quality 
of the exchange-correlation potentials. The first correction has been the 
Generalized Gradient Approximation (GGA), where a V^c is a function of 
the position, of the electronic density and of the density gradient. A proper 
GGA formulation, which respects the same internal sum-rules that are the 
key of the LDA effectiveness, was given by Perdew, Burke and Ernzerhof ''^ 
(GGA-PBE). In the so called meta-GGA exchange-correlation functionals 
also the density laplacian is included and in hybrid functionals an exact 
Hartree exchange term is considered. The quality of the improvements is 
not systematical, though: in some cases the LDA performs better than the 
GGA^; for this reason the research in this field is still very active. 

Of course, a DFT calculation involves a lot of technical details, related 
to the basis set, the pseudopotentials (when needed), the periodic conditions 
and so on. The specifics for the calculation of this Thesis work are given in 
the next paragraph. 
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2.6.3 Computing details 

The calculations of this Thesis have been done with the ABINIT pack- 
age®, based on a plane- wave basis set. Each KS single particle ip wave- 
function is thus described with a sum of plane- waves (PWs); this basis-set 
requires 3D periodic conditions and thus is best suitable for crystals, where 
the electronic states can be expressed with the Block theorem. The main ad- 
vantages of this method are the generality, (since the PW are unbiased) and 
the speed; in fact, with plane- waves the KS equations can be solved in the 
reciprocal space with the aid of the Fast Fourier Transform algorithm. The 
accuracy level of the results can be tuned by the energy cut-off parameter 
(Ecut) that fixes the total number of PWs. The complexity of the calculation 
scales with the real-space cell size. Molecules and surfaces can be simulated 
within larger supercells with a suitable region of empty space at the border. 

However, because of the diverging shape of the Coulomb potential, the 
core-level wavefunctions are highly oscillating around the nucleus and thus 
require a very large number of Fourier components to be properly described. 
For this reason, in a PW-DFT calculation the real ionic potential is usu- 



ally substituted with the so called pseudopotentials (see Fig. 2.6.3). A 
Pseudopotential (PSP) is specifically designed in order to produce wavefunc- 
tions that match the exact ones outside a certain distance from the nucleus, 
while a smooth non-physical solution is left inside. The pseudopotentials are 
non-local, since are different for each / quantum number. The diffusion of 
DFT calculation is related to the computational speed-up obtained by the 
introduction of PSP. The chemical properties are more dependent to valence 
electrons than core ones; therefore, a good description of the properties of 
molecules and crystals can be achieved by discarding a defined set of core 
electronic levels and by working with the few remaining electrons in a "pseu- 
doatom" with a lower ionic charge. 

The other possibility in the basis-set choice is to work with a localized 
set of functions; the most common are Slater-type Orbitals described with a 
linear combination of Gaussian peaks (STO-NG basis, where N is the num- 
ber of gaussian peaks). In transition metals it is possible to add a specific 
additional basis subset aimed to mimic the d orbitals. This basis set has the 
advantage of being localized around an atom and to have a relatively small 
number of parameters. Moreover, because of the gaussian peak properties, 
many of the electronic integrals can be calculated analytically. The disad- 
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Figure 2.15: Schematic view of 
a ionic (solid line) and a pseudized 
(dasiied line) potentials and wave 
functions, designed to match at a 
radius r^. Adapted from Ref. 'SD. 




vantage is the relatively fixed shape of the functions, which in some cases 
can poorly describe the long-range electronic states. The typical example 
are the carbon-nanotubes with small diameter, where the conductive states 
in the tube axis cannot be predicted by means of DFT with a localized ba- 
sis set'^^'. The STO-NG and the numerical equivalent (see for instance the 
SIESTA code^) basis sets are thus often used in quantum chemistry for 
molecules, since they don't require periodic conditions. 

Another opportunity is to mix both worlds by separating the basis set 
in two part: a localized (for core-levels) function set inside a critical radius 
(r^), and plane waves in the remaining space in order to describe the valence 
electrons. This field includes the muffin-tin (LMTO) based approach, the 
so called all-electron calculations (LAPW methods) and, to lesser extent, 
the Projector Augmented- Wave (PAW) method used in this Thesis work. In 
particular, in the PAW method a linear transformation that connects the 
all-electron wave function \E'n with a "soft" one \E'„ is defined through a set 
of projectors pf. 



I^n) = |^~n) + J2^<P^ ' 4>^) {p^\'^ n) (2.23) 

i 

where the all-electron 0j and the pseudized (pi partial waves are equal outside 
the PAW sphere with radius ic- The wavefunctions of Eq. 2.23| are given by 



the superposition of planewaves and a localized set of spherical harmonics 
functions. The PAW method is usually faster than the norm-conserving pseu- 
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dopotential method, because the loss of speed due to the higher complexity 
of the basis set is counterbalanced by a much lower Ecut- In this Thesis, 
the PAW formalism has been used because it also allows the introduction of 
localized Hubbard terms in specific orbitals. This is the formulation of the 
total energy in the LDA+U approach®: 

ELDA+uHr)] = EldaHt)] + EHub[n'^] - ^dc[n'l (2.24) 

where Efjub is the Hubbard correction added to the specific a orbital of atom 
/ and is an additional term that removes a double counting error on 
orbitals. When applied to TM oxides, the DFT+U theory allows a correct 
description of band-gap and magnetic ordering. 

Calculations in this Thesis have been performed on a quad-core desktop 
machine. Both norm- conserving pseudopotential and PAW basis have been 
used. For most of the elements some shallow core levels have been included 
in the basis set, in order to obtain results consistent with the experiments; 
for example, titanium ions have been described in a 3s^ 3p^ 4s^ 4d^ configu- 
ration, instead of considering only the 4s and 4d electrons. The convergence 
on total energy and k-grid density has been checked carefully, and every crys- 
tal cell has been relaxed to an equilibrium position. The Hubbard U values 
have been found in the literature; these are not completely empirical pa- 
rameters, since they can be evaluated theoretically with other DFT ab-initio 
calculations ® . 
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Chapter 3 

Ti02_(5 Resonant photoemission 



3.1 Introduction 

Titanium dioxide (Ti02) is one of the most studied metal-oxides. The 
constant interest on the Ti02 surface is due to the apphcative properties 
(especially as a single-crystal surface), to its simplicity in growth and in its 
electronic structure. The two most common Ti02 minerals are rutile and 
anatase, the former being the most abundant (stable) form. The rutile single 
cell is composed of six atoms (two titanium and four oxygen atoms) and 
belong to the tetragonal space group {Pi^/mnm, n°136); each titanium atom 
is at the center of a distorted octahedral cage of oxygen, with a D^h symmetry 



(see Fig. 3.1 ). 



The valence electrons are localized in the oxygen 2p levels, while the 
empty states are composed by Ti d levels. Ti02 should be a (formally simple) 
closed-shell d° band insulator; however, the hybridization of titanium 3d and 
oxygen 2p states can induce a partial occupation of Ti d-level. In fact, Ti02 is 
a charge-transfer insulator, like other TM oxides'^. The (direct) band-gap 
of rutile and anatase is of nearly 3.2 eV. 

The main applications of Ti02, such as photocatalysis, photodegradation 
of carbon compounds, solar cells and high-k gate for transistors®, require 
a fine tuning of the electronic structure, especially of titanium 3d states. 
In particular, the Ti02 based photocatalytic process require to shine the 
surfaces with UV light, in order to pump electrons from the top of the valence 
band to the empty states. A reduction of the Ti02 band-gap, achieved by 
selective doping, could in principle extend the catalytic properties up to the 
visible (solar) light range, with an obvious impact on the applications. 
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The understanding of the relationship between the defects electronic struc- 
tures and the so called "d° mae; netism" ® in TiOo poses an even more chal- 
lenging task, from both the theoretical and the experimental point of view. 
A detailed analysis of the magnetism in RF-sputtered anatase and rutile thin 
films is given in the next Chapter. 

The Resonant photoemission spectroscopy (ResPES) on transition met- 
als (TM) provides the opportunity to emphasize elemental-specific electronic 
structures in the gap and in the valence band. Theoretical prediction of res- 
onant behavior at 2p edge of early TM oxides have been done by Tanaka and 
Jo®, followed by the experimental work of Prince et al.''^ on Ti02 rutile. 
Resonant photoemission on different Ti02 phases have already been mea- 
sured on 3p-3d (or M3^4) edge, which gives higher sensitivity to surface 
states than the 2p-3d (L2,3) edge. While highly defective Ti02 films have 
already been reported®, data on non-conductive, slightly defected speci- 
mens are still missing, although recent studies on Ti-based oxide have shown 
the importance of the in-gap defect states in slightly substoichiometric sys- 
tems®. In this chapter, a ResPES study of transparent, non-conducting 
rutile Ti02-5 fil^i is presented. 



3.2 Resonant Photoemission on Ti02 



47 



3.2 Resonant Photoemission on Ti02 

As already explained in Chap. 2, in the case of TM compounds, the 
L-edges ResPES effects are expected to be rather pronounced. The Fano 
factor in TM-based systems is usually large (g > 3) and it will be shown, is 
even larger in Ti02. In our case study, photon energies have been scanned 
through the L2,3 edge; in a configuration-interaction approach the initial, the 
final and intermediate states in a general TM can be labeled as follow: 

2/3ci" -> 2p^3ci"-^ + e" (3.1) 



2p^3rf" ^ 2p^3ci"+^ (3.2) 



2p^3r^^ 2/3rf"-^ + e" (3.3) 



Eq. |3.1| defines the photoemission process, Eq. 3.3 defines the autoioniza- 



tion channel and Eq. |3.2| shows the 2p-3d transition. In the empty shell 
(3d°) Ti02 film, the VB resonances shouldn't be observed, because at least 
two d-electron are needed in the intermediate state in order to yeld the au- 
toionization process. However, in this case ResPES occurs because of the 
aforementioned hybridization of Ti 3d level with valence-band O 2p level. 
As shown in Chapter 2, this issue can be addressed with the charge-transfer 
multiplet theory and it is partly considered in DFT calculation (which already 
include some electron correlation, as defined by the exchange-correlation po- 
tential). 

3.3 Theoretical calculations for ResPES 

In order to theoretically calculate a ResPES spectrum, at least two differ- 
ent models have to be taken into account: an atomic charge-transfer model 
for the absorption cross section and an ab-initio density of states calculation, 
usually ground-state DFT. The first model generally includes several param- 
eters (charge-transfer energies, Hubbard U d-d repulsion, inter-configuration 



48 



Ti02 <5 Resonant photoemission 



mixing, crystal field) but can successfully describe matrix elements for ex- 
cited states, which determine the absorption edge on TM L-edges; the latter 
(DFT) doesn't need many empirical parameters (except U for DFT+U, which 
can be calculated from similar ab-initio approaches®') but it is difficult to 
extend to excited states. Usually''^, in a configuration- interaction approach, 
ground (g.s.) and final (f.s.) states of Ti02 resonant photoemission can be 
described as: 

= ao|3d°) + ai\?,d^l}) + aslSrf^L^) (3.4) 



■^f.s = (3o\2p3d^) + Po\2pM^L^) + Pq\2pU'^L') (3.5) 

where L and p stands for a hole in the ligand (O) and in the 2p (Ti) level. In 
rutile Ti02, an average occupation of d-levels of 0.7 can be estimated from 
the a coefficient, using literature data (q!o=39.5%, q;i=48.2%, Q!2=12.3%). 
The charge-transfer multiplet model has proved to be really effective for cal- 
culating XAS spectra at L2,3 edges of TM, because of the possibility to evalu- 
ate the transition matrix elements to the excited states; however, atomic ap- 
proaches fail in modelling VB spectra®, which must often be calculated us- 
ing ab-initio ground-state calculation, like density functional theory (DFT). 
An highly desiderable complete ab-initio approach is still a matter of re- 
search; good results have been achieved in modeling Ti L-edge XAS using 
multiple scattering theory®. 

In order to reproduce the non-resonant valence band spectra of Ti02, a 
DFT calculation have been carried out, both for a single rutile cell and a 



2x2x2 48 atoms supercell (with an oxygen vacancy, shown in Fig. 3.2); the 
calculations have been performed in the LSDA+U approximation with the 
ABINIT® code, using a projector augmented waves (PAW) basis set''^. 
The initial spin has been chosen to be up (down) for the first (second) Ti 
atom in the single rutile cell. The Hubbard U correction has been necessary 
to recover the correct band-gap (3.2 eV) and to obtain an insulating ground- 
state. In fact, an LDA (without U) calculation on this defected supercell 
usually gives a metallic behavior, as already reported in literature®. 

The LSDA-I-U calculation shows a net difference in spin-up/spin-down 
electronic densities, located at the three Ti atoms surrounding the oxygen 



3.3 Theoretical calculations for ResPES 



49 




Figure 3.2: LDA+U DFT Relaxed rutile 2x2x2 supercell, cutted on the plane 
containing the oxygen vacancy (titanium atoms are depicted in cyan, oxygen atoms 
in red); superimposed to the structure, a slice of the spin- up/spin-down difference 
electronic density is shown. In this calculation, U=8.0 eV (from Ref. ''^). 



vacancies. Every oxygen atom is bonded to three titanium atoms: the two 
Ti atoms belonging to the edge-sharing octahedra take a relative "plus" sign 
(yellow/hot colors in Fig. 3.2) while the apical octahedra has a "minus" 
sign (blue/cold colors in Fig. 3.2). These results, that concern the magnetic 
properties of Ti02, will be extensively discussed discussed in Chapter 3. 



From a technical point of view, energy cut off was 25 Hartree, and full cell 
relaxations have been performed with a maximum residual forces on atoms 
at 10~^ Hartree/A.; results are consistent with other all-electron calcula- 
tions''^. The X-ray absorption spectrum and core-level photoemission on 
Ti edge have been calculated with the MISSING package ''^3 by considering full 
multiplet and crystal field effects, in symmetry; crystal field parameters 
for this structure have been set at 10Dq=1.7 eV, ds=dt=-0.12 eV. 
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3.4 Experimental details and results 

The sample has been grown with RF sputtering technique from a Ti02 
target. The film thickness was thin enough (10 nm) to avoid charging effect. 
Resonant photoemission and x-ray absorption data have been collected at the 
BACH beamline of the ELETTRA synchrotron radiation source in Trieste 
(Italy). Usually, Ti02 surfaces for photoemission experiment are prepared by 
an heavy annealing in vacuum (600°C or higher), followed by an annealing 
in oxygen atmosphere to recover surface oxidation; this procedure usually 
generates a very ordered surface reconstruction and also make the Ti02 film 
conductive, allowing synchrotron measurement. Rutile crystals, after this 
procedure, usually change the color from transparent to black/blue, due to 
the annealing-induced defects. Since this work aims to measure the pristine 
crystal defect electronic structures, the Ti02 film has been left untreated 
(transparent and insulating); a layer of atmospheric contamination is thus 
expected on the sample surface. Even if this hydrocarbon contamination 
layer can attenuate the photoemission intensity, it can help to fix the binding 
energy scale (through XPS on C Is peak) and to saturate the dangling bond 
at Ti02 surface, so that only "bulk" Ti^~'^+ can be detected. 

The resolution of the monochromator was set to 0.24 eV at the Ti L2,3- 
edge photon energy. X-ray absorption spectroscopy (XAS) and photoemis- 
sion measurements have been carried out by using a modified 150 mm VSW 
hemispherical electron analyzer with a 16-channel detector. The total reso- 
lution for ResPES was 0.3 eV. XAS spectra have been taken in total yield 
mode, with circularly polarized light perpendicular to sample. During pho- 
toemission measurement, the sample has been kept normal to the analyzer. 



In Fig. 3.3 the measured and calculated XAS at Ti L-edge are shown. 



along with the Ti 2p XPS core levels on the thin film (Fig. 3.3, lower panel). 
Ti 2p core levels are compared to the results obtained from a high-quality 
single crystal rutile (110) surface. The XAS data compare well to the multi- 
plet calculations and are in good agreement with the results already reported 
in the literature for high quality samples®. The XAS spectrum shows a se- 
quence of four main bands, with alternating t2g and e^ symmetries (labeled in 
the figure) . The first couple of bands at low photon energies is separated from 
the second at higher photon energies by the spin-orbit interaction of the Ti 
2p core hole. The eg band at about 459 eV is actually split into two features, 
due to the crystal field effect that reduces the octahedral to a tetragonal sym- 
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Figure 3.3: Upper panel: Ti 2p-3d X-ray absorption spectra, with multiplet cal- 
culation in D4/J symmetry; lower panel: Ti 2p XPS spectra from the polycristalline 
film and from a reconstructed (100) rutile single-crystal surface. 
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metry; therefore, the D4/1 symmetry must be taken into account in order to 
reproduce the experimental data. The XAS data apparently do not suggest 
the presence of oxygen vacancies or Ti'^+ ions, which, in turn, can be related 
to the shoulder of Ti 2p3/2 XPS peak (black arrow in Fig. 3.3, lower panel). 
However, this finding does not allow us to draw further conclusion on the 
role of oxygen vacancies on the electronic properties of insulating rutile and 
therefore we resorted to ResPES techniques to enhance the spectral weight 
of Ti-related states in the valence band region, and to probe the relationship 
of these states with oxygen vacancies. 



In Fig. 3.4 the ResPES data on the Ti L-edge are shown. Each line of 
panel (a) represents a single photoemission spectrum taken at the correspond- 
ing photon energy, labeled on the vertical axis. The XAS profile across the 



Ti L-edge is also shown on the left axis. Fig. 3.4 c) shows the ResPES data 



obtained from panel (a) after the subtraction of the off-resonance spectrum 
collected at hu = 450.0 eV; this procedure has been done in order to en- 
hance the contribution of resonating spectral weight over the non-resonating 
valence band states. 

Selected valence band spectra are shown in Fig. 3.4[ d), while in Fig. 
|3.4K b) and in the inset of panel (d), enlarged views of the in-gap states, with 
BE ranging from to 2.5 eV, are shown. For a description of the overall fea- 
tures of the valence band photoemission region, we refer to Fig. 3.4[ d). Peak 
A is assigned to the oxygen 2s shallow core level. The valence band shows 
two features, D and D', the former resonating at hi'=458.1 eV. Therefore, 
the high energy states C clearly resonates at /iz/=458.1 eV, along with the 
in-gap state E, that resonates also at /ii/=456.3 eV. A weaker resonance is 
detected for peak B. 

The analysis of the DFT calculations for ground state of the 2x2x2 su- 
percell provides useful indications for the identification of some of the above- 
mentioned features. As it can be seen in Fig. |3.4[d), the LDA+U calculated 



density of states (DOS) can reproduce the shape of the valence band (VB), 
dominated by the oxygen 2p band, and provide information about the Ti-3d 
projected DOS in the valence band that, because of the O 2p - Ti 3d hy- 
bridization, results to be non negligible (especially in the high BE part of 
the VB, i.e. below the experimental D band). Therefore, peak D (BE=8.2 
eV) is the resonating structure related to the Ti-0 bonding part of the VB, 
while the lower BE region corresponding the peak D is usually referred to as 
the non-bonding part of the VB ^ . 
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Binding Energy (eV) 



Figure 3.4: RePES data: (a) contour plot of the raw data, along with the Ti 
2p-3d absorption spectrum (left axis); (b) enlarged view of the defect electronic 
states states nearby the Fermi level; (c) ResPES data after the subtraction of the 
off-resonance (/iz/=450 eV) spectrum; (d) selected photoemission spectra spectra 
(collected at /ii^=450, 456.3 and 458.1 eV) compared to LDA+U calculations of 
the total and Ti 3d projected DOS. Inset: details of the in-gap electronic states 
measured at /iz^=450, 456.3 and 458.1 eV. Taken from Ref. 'I^. 
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The low intensity peak (E) found at BE=1.5 eV, can be seen in detail 
in Fig. 3.4[ b) and in the inset of Fig. 3.4[ d) for selected photon energies. 
This peak, that is not detectable off-resonance, appears to resonate at differ- 
ent photon energies with respect to main resonance (maximum at hi'=A56.3 
eV rather than /it/=458.1 eV) and it is ascribed to oxygen vacancy states, 
yielding formally Ti^+ ions nearby the vacancy itself. Indeed, the DOS cal- 
culation on the defected supercell (Fig. 3.4[ d)) well reproduces the position 
of these in-gap states. Notably, the position of E peak seems to slightly shift 
upon changing photon energy from the pre-edge (/iz/=456.3 eV) to the main 
absorption edge (/ii/=458.1 eV), revealing a not completely resolved manifold 
of states. In fact, the slightly asymmetrical lineshape of E peak is similar to 
the one observed in other Ti-based 3d^ highly correlated systems (see, e.g., 
Ref. and Refs. therein). 

We now address the origin of peak B and C by referring to the available 
results of configuration interaction calculations of the photoemission process 
as well to the tracking of B and C peaks BEs upon photon energy increase. 
The |3c?"'Z/™)-like configurations of Eq. 3.4 are those accounting for charge 
transfer (CT) effects from the anion 2p to the TM 3d orbitals, providing 
a straightforward tool to consider the effect of 3d-2p hybridization in CI 
calculations. Therefore, the possibility to observe ResPES effects in the 
formally empty 3d-shell of the Ti^"*" cation is due to the presence of CT-like 
configurations both in the initial and final states. Furthermore, the presence 
of CT-like configurations in the initial state should agree with the finding that 



the Ti 3d projected DOS (Fig. |3.4[ d)) in the valence band is not negligible. 

As a first step of our analysis, we focus on the most intense peak C of Fig. 
|3.4K c). As it can be observed, the BE of peak C is constant for low photon 
energies, but it shifts to higher values as the photon energy increases above 
458 eV. Therefore, peak C shows a radiationless Raman-to-normal Auger 
transition''^ and above /ii/=458 eV, peak C ends up to be the normal LVV 
Auger. Hence peak C is assigned to \3(fL^) initial state configuration, as at 
least two d electrons are required to yield a normal L3VV Auger emission. 
Apparently, the dispersion on the BE scale of peak C is not strictly linear, 
since a sequence of absorption thresholds is accessed in a relatively narrow 
energy range, introducing new resonances that make the BE tracking of this 
peak rather difficult. However, it is also clear that above about /ii/=463 eV, a 
new normal- Auger feature (L2VV) becomes detectable, as the photon energy 
crosses the Ti L2 threshold. Finally, peak B also shows a dispersion above 
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the resonance on the BE scale, though weaker than that observed for peak 
C. 

It is important to note that, unlike peak C and, to a minor extent, peak 
B, peak D does not show any detectable dispersion on the BE scale. There- 
fore, in the CI approach, for peak D a major contribution of the |3(i^L^). 
configuration is expected, as for this configuration the autoionization chan- 
nel is allowed, but not the normal Auger emission. According to Runs et 
al.''^, peak D is fit by the most intense feature of the CI calculations®, 
and starting from this peak all the higher-BE peaks (B and C) are well fit 
by the CI model. On the other hand, the energy of peak D corresponds to 
the region where the Ti-projected DOS has the largest intensity, therefore 
providing an alignment for the two computational results. However, while 
the single-configuration, ground-state, LDA+U calculations fail to reproduce 
the high BE satellites B and C, this indicates that the Ti 3d states are not 
only localized below the peak D, but they are spread, through band structure 
effects, to lower BEs, where CI calculations do not predict any Ti contribu- 
tion, as the band dispersion of oxygen states is usually not included in the 
CI model. 

In Fig. |3.5K a)-(c) an enlarged view of the ResPES data collected in the 
low BE range (0-10 eV) is shown, while Fig. 3.5[ d) shows the constant 
initial state (CIS) profiles of peak D, (BE=8.2 eV), peak D' (BE=6.0 eV) and 
defect states (peak E at BE=1.5 eV). These intensity profiles are obtained by 
vertically slicing Fig. 3.5 |^a), along the three dashed lines at constant BE. Due 
to a very high Fano q-factor (q>10), at the TM L-edge, the CIS profiles of 
selected valence band states usually closely follow the shape of the absorption 
edge. However, the CIS spectrum for the defect structure (BE=1.5 eV) is 
clearly different from the Ti L-edge XAS spectrum, especially in the case of 
peaks A, C, and F of Fig. 3.5[ c). This CIS spectrum could be regarded as 
an effective on-site defect XAS and it is compatible to that found in heavily 
sputtered rutile Ti L-edge XAS data®. The narrow width of the present 
CIS peaks with respect to those observed in heavily defective samples (see 
e.g.^) could be related to the reduced lattice disorder effects with respect 
to the heavily sputtered surfaces. 

Finally, the ResPES spectrum collected with a hu = 461.5 eV photon 
energy shows a clear resonant enhancement of peak D', that seems to shift 
to higher BEs while increasing the photon energy. In fact, a careful tracking 
of peaks BE reveals the presence of three Auger-like features that start to 
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Figure 3.5: Detail on valence band part of ResPES: in panel (a), ResPES on 
VB without the out-of-resonance contribution;in panel (b), spectra taken from 
image of panel (a) together with DFT m-resolved d-DOS calculations; in panel 
(c), ResPES normalized to the absorption edge. 
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resonate in the non-bonding part of VB at plioton energies 461.5 eV, 463.5 
eV and 466.2 eV (dotted lines in Fig. 3.5[ a)). 

To further explore this effect, in Fig, 3.5[ b) the VB ResPES data normal- 
ized to XAS spectrum intensity are shown: by doing so we can quench the 
XAS matrix contribution and enhance the spectral weight due to interfer- 
ence effects in the ResPES process. Three spectral features, Al, A2 and A3, 
displaying a BE shift appear, Al being the most intense (marked by arrows). 
These features are too prominent to be assigned to pure 3d^ defect states, 
and resonate mainly around the e^ symmetry edges of the absorption, while 
the more intense LVV Auger peaks are generated on t2g edges. Therefore, 
this state (peak D') presents a normal Auger (i.e. a two hole correlated 
satellite) behavior similar to the one observed for peak C, indicating that in 
peak D' a large contribution of the \3d'^l/) configuration in the initial state 
is expected. Incidentally, we recall that the e^ XAS bands are those mostly 
affected by the CP effects arising from the D^h distortion of the TiOe clus- 
ter cage. The symmetry dependence, which is ultimately a dependence on 
the orbital degrees of freedom, can also be discussed by considering the pro- 
jected DOS in the VB region as in Pig. |3.5K b), where the m-resolved d-DOS 
(m=0,±l,±2) are shown. While the strongest contribution to the bonding 
part of VB comes from the the m=0 and m=±l, related to d^^, dxz and dyz 
orbitals, the m=±2 projected DOS, which is related to dxy, dx2-y2 orbitals, 
provides a relevant contribution to the energy region below peak D'. 

In conclusion, a ResPES experiment at the Ti L-edge has allowed us to 
probe the in-gap states of an insulating Ti02-5 thin film, otherwise missed 
by off-resonance measurements. The corresponding CIS profile indicates that 
these states should be ascribed to Ti^+ ions. In spite of the rather complex 
sequence of thresholds, determined both by crystal-field and spin-orbit inter- 
actions, by tracking the dependance of the BE of each resonating peak on 
photon energy we have been able to identify the main configurations con- 
tributing to these electronic states. Distinct two-hole correlated satellites in 
the final state are related to the spectral features at BE=6, 13 and 19 eV. 
In particular, we provide evidence of a character for the peak below 

the overwhelming non-bonding states of the VB (D' at BE=6 eV), so far 
virtually neglected by CI calculations, as well as by the analysis of exper- 
imental data. Our finding asks for more refined CI calculations aimed to 
draw a consistent picture of the N-1 electron final state of the photoemission 
process. It is likely that this split-off peak could be properly described by 
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(i) explicitly accounting for the oxygen 2p band width in CI calculations or 

(ii) by resorting to a TiOe cluster symmetry lower than the one (Oh) so far 
assumed in CI calculations for the photoemission process I^SIMJ ^ agreement 
with the symmetry (D4/1) already adopted for the calculation of the Ti L-edge 
XAS spectrum. 



Chapter 4 



Magnetism in Ti02_5 
N:Ti02_^ films 



4.1 Introduction 

The history of diluted magnetic oxides (DM0) started in 2000 with a 
work of Dietl''^, in which the authors theoretically predicted an high-Tc 



ferromagnetism (FM) in Co-doped ZnO (see Fig. 4.1). Such kind of mate- 
rial would display a simultaneous presence of magnetism and optical trans- 
parency, due to the high band-gap, and would be a very interesting system 
for spintronic applications. Subsequent experimental work didn't find FM in 
ZnO-based compound*^, but many research group reported the presence of 
room-temperature FM in Co-doped Ti02 film ''^2'^. The Curie temperature 
in this cases is expected to be higher than 500 K; however, this kind of FM 
is usually rather weak (in term of coercive field) and thus it is not suitable 
for typical applications. 

The growth process has proved to be crucial for the magnetism in DM0. 
For instance, PLD (Pulsed Laser Deposition) grown Ti1_a.C02.O2 samples dis- 
plg^ydUli ^ very high magnetic moment per Co atom (1.5-1.7 ^b) which is not 
expected in a diluted sample and is caused by metallic cluster precipitation. 
Films of the same material, grown by reactive sputtering '^3' , display a much 
lower value (0.96 (Ib), that could be consistent with a non-clustered random 
Co distribution or with a mixture of metallic cluster and isolated Co atoms. 

The difference between an "intrinsic" or "extrinsic" origin of the mag- 
netism is fundamental for the theoretical models and for real applications. 
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Figure 4.1: Theoretical prediction of ferromagnetic Tc in different materials, 
upon doping with Mn ions. The wide-gap material are predicted to give higher 
Tc. Taken from Ref.lSS. 



Only an intrinsec magnetism would be useful for the spintronc field, since the 
charge carriers have to lie in a spin-polarized band. It is generally believed 
that it is possible to grow diluted Tii_^. 00^^02, with a suitable control of the 
growth parameters; unfortunately, in the attempt to achieve the highest Tc, 
many research groups focussed their work only on cobalt doping, which can 
easily generate secondary magnetic phases (and thus extrinsic FM). 

Subsequently, FM has been reported also in sing le crystal Gi» TM-doped 
Ti02 (TM = V, Cr, Mn, Fe, Co, Ni, Cu); the magnetization measurements 
appear to be similar for each dopant, with a different saturation magneti- 
zation. Other research group also reported FM on TM-Sn02 and TM- 
ZnO® systems. Therefore the driving mechanism of such kind of FM, 
observed for every kind of TM doping (some of which usually don't develop 
FM phases with oxygen, like Mn), has to be related to a general mechanism 
instead of a specific dopant. The magnetism is then likely to be mediated by 
intrinsic defects, among which the most common in these systems are oxygen 
vacancies (Vo's). 

In fact, one side effect of the doping is the replacing of Ti^+ ions with 
stoichiometric non-equivalent ions that can induce structural defects. In this 
description, the presence of magnetism should be triggered by interacting 
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Ti^"*" ions, induced by Vo's; tlie dopant role should be both to create Vo and 
to increase the saturation moment. 

These hypothesis has been further confirmed by the finding of ferromag- 
netism in undoped Ti02-5®, Hf02''^ and other® closed-shell oxides (CaO, 
ZnO, SrTiOa, MgO, LaAlOg, LSAT, AI2O3) and non oxide materials (AIN, 
GaN, proton-irradiated graphite); these materials are classified as "d° mag- 
nets". In each of these compounds, a sort of structural disorder is needed 
to observe the FM '^3 and the required density of defects is usually achieved 
with non-equilibrium growth methods. Moreover, the measured FM in d*^ 
thin films is usually highly anisotropic, being much stronger with an exter- 
nal field parallel than perpendicular to the surface. 

From a theoretical point of view, DM0 's and d° FM poses several prob- 
lems; first of all, magnetic interaction appears at doping-level far beyond 
the percolation threshold in the lattice: in other words, assuming an homo- 
geneous distribution of magnetic centers, dopant ions (or Vo) should not 
interact, because their average distance is far too high for super-exchange 
interaction to take place. The RKKY magnetic coupling mechanism is ex- 
pected to work only in the presence of a metallic conduction. A "localized" 
mediator to magnetic interaction could be given by trapped electronic states 
at Vo sites, while cationic vacancies are known to carry a finite magnetic 
moment due to the open shell structure'-^-'. The kind of defect that is re- 
sponsible for the magnetism is still unknown. 

From an experimental point of view, dealing with crystalline defects in 
oxides can be rather challenging. The usual way'^ to characterize the de- 
fects in Ti02 is to perform measurements (in particular, transport, XPS and 
EPR) on highly annealed samples (up to 1300 °C) that cause the oxygen 
desorption. Upon annealing, transparent insulating Ti02 crystals usually 
undergo a pronounced color change (up to dark blue) and to an increased 
conductivity. Annealed and re-oxidized films are often used in Ti02 photoe- 
mission surface studies in order to avoid charging effects, assuming that the 
annealing in oxygen could restore the full stoichiometry at the surface (for a 
thickness larger than the probing depth of the selected technique). However, 
it has been pointed out that bulk defects can play an important role in the 
surface structure® and, for example, that the sample preparation is thus 
crucial when performing resonant photoemission experiments''^. 

The presence of oxygen vacancies, interstitial Ti, clustered defects and 
planar dislocations have been observed''^, but the question of which kind of 
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Figure 4.2: The TiOx phase diagram as a function of x and of the anneahng 
temperature. The vertical hues in the 1.5 < x < 2.0 range correspond to the 
Magne h phases (Ti„02n-i)- Taken from Ref. ' — 



defect is prevalent for a defined stoichiometry is still under debate. For ex- 
ample, it has been proven that interstitials Ti (and not the oxygen vacancies) 
migrate on the surface during the re-oxidization process*^. The Ti-0 phase 
diagram, as a function of the Ti-0 ratio, is rich of secondary phases (see Fig. 



4.2); in particular, some elements of the homologous Ti„02n-i series of the 



Magneli phases (indicated by the vertical lines on the right part of Fig. 4.2), 
ranging from Ti02 to Ti203, can display magnetic properties. 

Furthermore, these magnetic effects are sometimes so weak (in terms of 
saturation magnetization) that a lot of effort has to be made in order to 
control the possible external contamination on the samples properly. For 
instance, a clean ceramic sample heater has to be used during the growth 
and post-annealing treatment, since a metallic one could cause metallic ion to 
diffuse into the sample. The samples have to be handled with (clean) ceramic 
or titanium tweezers, and can't be mounted on sample-holder with adhesive 
tape in order to avoid external contaminations. The XRD spectra have to 
be collected carefully in order to evidence even small peaks due to metallic 
inclusions. When working with powders or deposition targets purchased by 
chemical industries, an evaluation of the impurities at the ppm range should 
carried out, in order to cross-checking the declared composition. 

It is generally believed that, taking the necessary caution, a ferromagnetic 
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phase is really present in these oxides. Even if the possible application of 
these materials seems unlikely (despite the high Tc) the puzzle of this kind of 
magnetism is challenging from both the experimental and theoretical point 
of view and deserves further investigations. 

In this chapter, a series of magnetic measurements on several Ti02-5 
thin films is given. Hysteresis loops have been recorded with SQUID on 
"as deposited" sample, after an annealing in vacuum (in order to induce 
more Vo), and finally after an annealing in a O2 atmosphere (in order to 
recover the stoichiometry) . Samples have then been characterized with XPS 
on Ti 2p peaks but, as will be shown in the next sections, the effect of such 
small quantities of Vo's can't be detected with normal photoemission and 
ResPES seems to be the proper technique to address this problem. Thus a 
comparison of ResPES data on rutile Ti02-5 and N-doped Ti02-<5 thin film 
is given. Experimental results are compared to DFT+U calculation on Ti02 
and N-Ti02 supercells in the last chapter's section. 



4.2 Growth method and characterization 

The titanium oxide films have been grown with RF-sputtering from a 
Ti02 target on quartz substrate and, in order to improve the sample crys- 
tallinity, on oriented LSAT(OOl) and Al2O3(1102) surfaces. On the quartz. 



the control over the Ti02 crystal phase (anatase or rutile, see Fig. 3.1) was 
achieved by checking the temperature and the Ar/02 ratio in the gas mix- 
ture. The growth on perovskite LSAT ((LaA103)o.3(Sr2AlTa06)o.7) results in 
fully Anatase phase samples and the growth on sapphire in the rutile phase, 
without the need of additional O2 gas during the growth. While rutile is the 
most thermodynamical stable crystal structure, often thin films are grown 
in the anatase form; both are expected to display ferromagnetism. In order 
to evidence interface or surface effects, a "thick" and a "thin" (less than 60 
nm) film have been grown for each substrate. A list of the samples is given 
in Table 14.11 

The samples grown in an higher Ar pressure (< 2-10~^) or in a lower 
pressure 650°C substrate temperature resulted in a mixture of rutile and 
anatase phases. The magnetization of each sample has been measured with 
a SQUID magnetometer after the growth process (labeled deposited" 
films), after a 4 hours annealing in vacuum (10~^ mBar) at 400 °C ("re- 
duced" films) and after the same annealing treatment in a O2 atmosphere 
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Substrate Temp. Pressure ©2% Phase Dep. Time Thickness 





°C 


mBar 






min 


nm 


Quartz 


650 


5-10-3 


9 


Anatase 


45 


27.4 


Quartz 


650 


3.7-10-3 


9 


Anatase 


165 


83 


Quartz 


750 


4-10-3 





Rutile 


30 


57.7 


Quartz 


750 


4-10-3 





Rutile 


120 


226 


LSAT 


750 


5-10-3 





Anatase 


135 


317 


LSAT 


750 


5-10-3 





Anatase 


20 


47 


AI2O3 


750 


5-10-3 





Rutile 


135 


220 


AI2O3 


750 


5-10-3 





Rutile 


20 


33 



Table 4.1: List of Ti02 samples grown with RF-sputtering for magnetic mea- 
surement. 



("oxidized" films). It has to be noted that the annealing temperature was too 
low to induce structural modifications: the deposited films of table |4.ll were 



perfectly transparent and insulating, without any observable color change. 

Since photocatalysis is the most important Ti02 application, a nitrogen- 
doped rutile Ti02 sample has also been grown. Following the suggestion of 
Elfimov et al. an increased magnetism could also be expected in oxides 
with substitutional nitrogen ions; a structural and magnetic characterization 
has then been performed also in this case. The details of the growth param- 
eter and the SQUID measurement are given in the dedicated section of this 
Chapter. 

The crystalline order has been checked with XRD using a Bruker D8 
Advance diffractometer and with Micro-Raman using a Labram-Dilor spec- 
trometer. The sample growth, the magnetic and the structural character- 
ization have been done at the Universita degli Studi di Pavia. A Raman 
study on several Ti02 samples grown in the same laboratory has already 
been published (see Ref. ''^). 



4.3 Magnetic characterization 

Static molar magnetization (Mmoi) loops were collected at room temper- 
ature for magnetic fields (H) ranging between and ±4000 Oe with a super- 
conducting quantum interference device (SQUID). Reference measurements 
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have been done also on the SQUID sample holder and on the substrates, 
in order to evaluate only the contribution of the deposited films. A typical 
plot of the magnetic moment versus the external magnetic field is shown in 
Fig. 4.3, which represents the raw data as measured by SQUID on the thin 
anatase sample on quartz. 
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Figure 4.3: Magnetic moment vs external field measured on the thin anatase 
film deposited on quartz, taken at T=300 K. This graph shows the "raw" SQUID 
data, with the ferromagnetic and paramagnetic contribution. 



A hysteresis cycle is detectable both in the as-deposited and treated 
films and it is superimposed to a paramagnetic (PM) (or diamagnetic) linear 
background. In fact, high-quality single-crystal Ti02 should display diamag- 
netism (DM), which is indicated by a negative slope of M(H); the presence 
of oxygen vacancies due to the growth or induced by the vacuum annealing 
may yield a paramagnetic contribution that reduces or even eliminates the 
diamagnetism. The hysteresis loop of each data-set was extracted by remov- 
ing the background obtained trough a linear fitting on the high magnetic 
field part of spectra. 

A summary of magnetic measurement of Ti02 film grown on quartz sub- 
strates is shown in Fig. 4^ The data are already separated in the ferro- 
magnetic and paramagnetic (or diamagnetic) contribution. The as-deposited 
thin rutile film is diamagnetic (as expected for Ti02), while the thin anatase 
and thick rutile display paramagnetism, which indicates the presence of de- 
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Sample 


^as—dep 


^reduced 


^oxidized 


Thin anatase 


0.037 


0.312 


0.074 


Thin rutile 


0.015 


0.059 


0.044 


Thick rutile 


0.052 


0.014 


0.007 



Table 4.2: Value of 6 in the Ti02-5 films grown on quartz, estimated from the 



value of the saturation magnetization, using the Eq. 4.1 



fects due to growth. The green arrows shown in Fig. 4.4 a)-(b) underline the 
effect of the first in-vacuum annealing, which increases both the saturation 
magnetization and the paramagnetic contribution, the latter identified by a 
counter-clockwise rotation of the linear background. In all cases, the subse- 
quent annealing in oxygen, marked by the red arrows, restores (or improves) 
the crystalline order of the all samples, quenching the ferromagnetism. 

The samples have been carefully manipulated with ceramic tweezers, and 
the annealing has been done with a proper ceramic heating stage. In fact, 
assuming that the magnetism were due to the diffusion of FM particles from 
the heater, the magnetism should have increased also after the second an- 
nealing (the one in oxygen atmosphere), which is not the case. The data 
shown in |4.4| thus clearly display the importance of oxygen vacancies in the 
magnetism of Ti02. 

It is possible to evaluate the number of magnetic ions per unit formula 
from the saturation magnetization. The total molar magnetic moment is 
given by the following expression: 

M = rifiegS (4.1) 



where n is the molar density of magnetic ions, g the Lande factor and S the 
spin. In the case of a Ti^+ ion, the expected spin is 1/2. In table 4.2 the 
evaluated film stoichiometry in term of Ti02-<5 is shown, extracted from the 



data of Fig. 4.4 



Apart from the reduced thin anatase (with a theoretical TiOi.7 compo- 
sition, in which a presence of Magneli secondary phases is expected), the 
values displayed in table 4^ are exceedingly low respect to the usual density 
threshold associated to magnetic properties in diluted samples. 

In fact, for a random distribution of dopants, a long-range order (like FM) 
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Figure 4.4: Magnetic moment vs external field measured on the Ti02 films 
deposited on the quartz substrate; in (a)-(b), the thin (27.4 nm of thickness) 
anatase sample; in (c)-(d), the thin (57.7 nm) rutile; in (e)-(f) the thick (226 
nm) rutile. On the left panel are shown the ferromagnetic contribution, on the 
right panel the paramagnetic or diamagnetic one. The effect of the first (second) 
annealing is marked by the green (red) arrows on panels (a)-(b). 
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appears when the size of the larger cluster of interacting atoms is comparable 
to the size of the material. At fixed interaction radius, the density for which 
the long-range interaction appears is given by the solution of the so-called 
percolative problem. The percolative threshold depends on the system ge- 
ometry (lattice type and dimensionality) and represents a way to simulate 
the effect of disorder in diluted system. The the magnetism of a Mn doped 
CdTe can be fairly understood by the percolation of interacting Mn ions in 
the CdTe lattice (see Ref.GS and Ref.^S for a 2D case). Values for the 
percolative threshold in various 3D lattice (including anatase) are given in 
the appendix of Ref . ^ . 

Indeed, the densities of magnetic centers listed in Table 4^ are much lower 
than the percolation threshold in rutile or anatase lattice. This is a common 
result of d° ferromagnetic material, and justifies the cautious behavior of the 
research community about claims of FM ordering. Surprisingly, the thick 
rutile sample seems to display (in average) a smaller ferromagnetism respect 
to the thin one. As will be shown, this behavior can be found at higher 
extent in the oriented substrate films. 

In principle, it should be possible to evaluate the number of magnetic 
centers that contribute to paramagnetism (or diamagnetism), by fitting the 
slope of the linear contribution with a Curie-type law; however, as paramag- 
netic effects can also be found in the substrates after annealing treatments, 
this evaluation is not reliable for the present systems. 

The SQUID data for the anatase films grown on LSAT, already separated 
from the linear background, are shown in Fig. 4.5, while the data for rutile 
films on alumina are shown in Fig. 4^ The results are comparable with the 
ones obtained on the samples grown on quartz. The annealing procedures 
results in the expected behavior: the in-vacuum annealing increased the FM 
by increasing the number of Vo's, while the oxidation resulted in the opposite 
effect. The absolute saturation magnetization and the relative stoichiometry 



in terms of Ti02-5 are listed in Table 4.3 



The extrapolated density of magnetic ions is of the same order of the one 
measured in the Ti02 grown on quartz substrate and thus is again too low 
to justify a long-range magnetic ordering. In particular, in an anatase- type 
lattice, for a cation-cation interaction the percolation occurs when 27.9% of 
Ti atoms have a magnetic moment, by considering only the nearest-neighbor 
interaction''^. According to Coey et al. ®, the high dielectric constant of 
these materials can increase the range of magnetic interaction between ions; 
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Figure 4.5: Magnetic moment vs external field measured on the thin(a) and 
tliick(b) film grown on LSAT. In black the as-deposited samples, in green the 
reduced one and in red the oxidized films. 
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Figure 4.6: Magnetic moment vs external field measured on the thin(a) and 
thick(b) rutile film grown on Al2O3(1102). In black the as-deposited samples, in 
green the reduced one and in red the oxidized films. 
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Sample 


^as—dep 


^reduced 


^oxidized 


Thin anatase 


0.031 


0.035 


0.003 


Thick anatase 


0.002 


0.003 


0.0003 


Thin rutile 


0.051 


0.136 


0.034 


Thick rutile 


0.078 


0.088 


0.040 



Table 4.3: Value of 6 in the Ti02-(5 films grown on LSAT and alumina, estimated 



from the value of the saturation magnetization, using the eq. 4.1 



in the case of Hf02, the cation-cation interaction can be effective up to the 
fifth-neighbor, which correspond to a doping concentrations of 13.7%, still 



higher than the values shown on Table 4.3 



Nevertheless, the absolute value of the saturation magnetization of our 
film becomes more significative when compared to a surface or to a 2D-like 
distribution of spins instead of a 3D one. A proof of a not-homogeneous de- 
fects distribution can found by looking at the saturation magnetization in the 
anatase thick sample of table 4.3 which is surprisingly low (0.25 emu/cm^) as 
compared to the thin one. A plot of the estimated number of Bohr magneton 



per unit formula versus the film thickness for each sample is given in Fig. 4.7 



The scattering of the data reflects the disordered nature of the phenomenon. 
In general, the thin films (especially after the annealing in vacuum) display 
a higher magnetization than the corresponding thick one. This should not 
happen in homogeneous diluted ferromagnet, in which the magnetization per 
unit mass (or volume) is constant. 

It should be pointed out that the magnetization has been calculated as- 
suming a uniform distribution of magnetic moments. The difference of m^at 
of the thin and the thick film can be partly reduced when the films thickness 



listed in table 4.1 is considered. An accumulation of defects (Vo, Vri or 
interstitial O and Ti ions) near the surface or the interface can justify the 
presence of FM interaction, since the higher local magnetic moment density 
could reach the percolation threshold. For example, according to the per- 
colation threshold (27.9%) estimated for nearest- neighbor interaction in the 
anatase lattice''^, in the reduced thick anatase film grown on quartz this 
threshold can be achieved if the total number of magnetic ions were confined 
in a 1.7 nm thick layer. Furthermore, the theoretical explanations in liter- 
ature (see Ref. ®), discussed at the end of this chapter, suggest that FM is 
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Figure 4.7: Number of Bohr mag- 
neton per unit formula versus the 
film thickness for each sample. The 
marker label correspond to the sam- 
ple: Q(0) is for sample grown on 
Quartz (oriented surface); A(R) is 
for Anatase (Rutile) film. 

caused by clustered defects. This model could also explain the anisotropic 
magnetism of d° thin films. 

The location of this "magnetic" layer should be near the boundaries of 
the film, i.e. at the surface or the interface. The reduction of the samples, 
induced by the annealing, occurs at the surface; however, the magnetism 
appears also in as-deposited samples, and the lattice mismatch and roughness 
could induce an higher quantity of structural defect at the beginning of the 
sputtering process, i.e. at the interface. Moreover, as pointed out in Ref. "I^, 
during the annealing process the reduction of Ti02-a: is accompanied by a 
migration of Ti atoms in the bulk that can increase the density of magnetic 
moment carriers at the interface. 

4.4 XPS measurements 

An XPS analysis has been carried out in order to investigate the pos- 
sible presence of defects at the surface. The thin and thick rutile sample 
grown on Quartz and on Sapphire substrate have been measured, both as- 
deposited, oxidized and reduced. The XPS measurements have been taken 
with a VG Escalab Mk II apparatus with a modified 32-channel detector 
and a monochromatized Al ka x-ray source in the Surface Science lab at 
the Universita Cattolica of Brescia. The charging effects induced by x-ray 
exposure have been compensated with an e~ flood-gun. Since the purpose 
of this study was the characterization of the defects electronic structure, we 
didn't perform further annealing and sputtering treatment. For transparent 
and insulating of Ti02 samples, i.e. for such a low amount of oxygen vacan- 
cies, it is not possible to directly detect the defect-induced modification in 
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Figure 4.8: XPS spectra of Ti 2p core levels on the thick Rutile film grown on 
quartz. 



the valence band. However, it should be possible to evaluate the presence of 
Ti'^"'" states by analyzing the photoemission structure of the Ti 2p core-levels. 
The results for the thick rutile sample grown on quartz is given in Fig. |4.8 



The three spectra of Fig. |4.8 are almost identical; it is not possible to 
clearly detect the presence of Ti^"*" related structures, which should appear 
at lower BE respect to the main peak (2p3/2, at BE=458.5 eV). The same 
results have been obtained for the other samples of Table |4.1 Assuming 
an homogeneous distribution of defects, with the densities listed in Table 
4]2] and |4.3[ a difference in spectra should be easily detectable. The typical 
probing depth of XPS at Ti 2p core-level, with a photon energy of hu =1486.6 
eV is given by the electron inelastic mean free path A =2.212 nm, obtained 
from the TPP-2M formula''^. Therefore, a difference in the Vo density at 
a distance A from the surface should be excluded in these Ti02 film. 

A more bulk-like technique that can prove the relative disorder of a crys- 

(Ref. 'l^) shows the peak 



4.9 



tal structure is the Raman spectroscopy. Fig. 
width of selected Raman peaks taken on films grown on quartz, as a func- 
tion of the sample treatment. The Raman modes of fully oxidized samples 
display a peak width that is considerably lower than the as-deposited or an- 
nealed films, as a consequence of the improved crystalline order. These data 
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Figure 4.9: Peak width of selected 
Raman mode for anatase and ru- 
tile sample grown on quartz. The 
in- vacuum annealed samples are la- 
beled as "vac", the O2 annealed as 
as-dep oxy vac "oxy" . Taken from Ref. . 

confirm the effectiveness of the anneahng procedure in the creation of de- 
fects, even if the anneahng temperature (400 °C) is considerably lower than 
the one typically used to induce oxygen vacancies in literature. Together 
with the XPS data, Raman spectroscopy could support the hypothesis of 
a localized layer of defects near the interface. The experimental noise and 
the difficulty of distinguish a small amount of oxygen vacancies in the XPS 
should also be taken into consideration. We thus resort to Resonant Photoe- 
mission (ResPES), in order to probe the electronic structure of even a small 
concentration of defects near the surface. 

Since the ResPES technique requires a synchrotron facility to be per- 
formed, and since a typical synchrotron X-ray beam causes much stronger 
charging effects than normal XPS, only the thinnest rutile Ti02 sample has 
been probed. The results of this analysis is given in the next sections and 
have been published in Ref. ®. 

4.5 ResPES on TiOs-^ and N-doped TiOs-^ 

As already pointed out, the typical densities of defects measured in fer- 
romagnetic Ti02 film is usually too low to explain a long-range magnetic in- 
teraction. This is particularly true for Ti02 where highly localized deep-gap 
defect orbitals lead to only a very-short-range defect-defect interaction. For 
very low densities of defects, these vacancies would fail to percolate through 
the sample, resulting in a vanishing Curie temperature. On this basis, a 
possible explanation for the origin of FM coupling in insulating DM0 can 
be the inhomogeneity of Vo distribution that can be high enough in limited 
regions to give FM coupling. The point is now to establish whether doping 
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Figure 4.10: Picture of the rutile Ti02 (a) and N-doped Ti02 (b) rutile film; in 
panel (c), the optical transmittivity; in panel (d), the optical reflectivity. 



of rutile Ti02 can create better conditions for the insurgence of FM ordering. 

In this section, we show the effect of N-doping on the magnetic and elec- 
tronic properties of Ti02-5 rutile thin films. The thin films (about 40 nm 
thick) have been grown by RF-sputtering a Ti02 target on oriented Alu- 
mina (0001) substrates at 750 °C under a flux of pure Ar (undoped film) 
or an ArN mixture (95%5%, N-doped film), at a pressure of 5x10^'^ mbar. 
The crystalline phase of each sample was checked by X-ray diffraction and 
micro- Raman spectroscopy. While the undoped film is perfectly transparent, 
the N-doped film display a yellow-gold color (see Fig. 4.10[ a)-(b)), that can 
be probed with simple optical transmittivity and reflectivity experiments, 
shown in Fig. 4.10[ c)-(d). 

The color change reflects the band-gap reduction, from nearly 3.2 eV 
to 3.0 eV, marked by the arrows in Fig. 4.10[ d). Many experimental re- 
search groups already reported an increased photocatalytic ''^ activity of 
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Ti02 under sunlight, achieved trough a band-gap reduction upon nitrogen 
or vanadium doping. 

Upon N-doping an enhancement of ferromagnetic ordering in terms of 
remanent and saturation magnetization is also achieved. By combining soft 
X-ray spectroscopy and magnetization measurements with ground state den- 
sity functional calculations, we are able to show that, with respect to the 
undoped system, N doping does not add electronic states in the region where 
Vo states are usually found, but introduces additional N-related low-lying 
acceptor and donor levels through O substitution with N. This ultimately 
reduces the band gap of rutile providing more favorable conditions for the 
onset of robust FM ordering in these compounds. 

Resonant photoemission spectroscopy (ResPES) and X-ray absorption 
spectroscopy (XAS) spectra have been collected at the BACH beamline of 
the Elettra synchrotron in Trieste (Italy). Ground-state density functional 
theory (DFT) calculations have been carried out, based on the Perdew- 
Burke-Ernzerhof (PBE®) parametrization of the generalized gradient ap- 
proximation (GGA) for the exchange-correlation functional. The projector 
augmented- wave ''^ atom description is used, as implemented in the ABINIT 
(Ref.®J) code, adding semicore level to valence states. Calculations have 
been carried out on a rutile cell (6 atoms/cell) and a 2x2x2 (48 atoms) su- 
percell with one nitrogen atom substitutional to oxygen, corresponding to a 
x=3.2% doping level. The k-space grid was a 4x4x4 MonkhorstPack grid, 
and the plane-wave energy cutoff was 25 hartree. We performed a full cell 
relaxation with a maximum residual force on atoms fixed at 10~^ hartree/A. 

Mmoi versus H curves collected at RT with the external field H applied 
parallel to the thin film surface (Fig. 4.11) show an hysteresis cycle that 
mark the presence of a ferromagnetic-like ordering already at RT. For the 
N-doped sample, saturation is reached at about 30 emu/cm^, a value much 
higher that the case of the undoped sample (about 5 emu/cm^). 

The effect of N on the electronic properties of the doped sample is dis- 
cussed by examining the Ti 2p and N Is core level x-ray photoelectron spec- 
troscopy spectra of the undoped and N-doped thin films. Figure |4.12[ left 
panel shows the Ti 2p XPS spectrum obtained from the undoped Ti02 thin 
film and from the N-doped thin film. These spectra are compared with the 
Ti 2p XPS core lines of a Ti02 rutile single crystal. The Ti 2p spectra show 
the typical spin-orbit split doublet of the Ti^"*" ion at a binding energy BE 
=458.8 eV and BE=464.8 eV. The Ti 2p3/2 line shows a shoulder at low 
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Figure 4.11: Hysteresis loop measured on the undoped Ti02 and on the N-doped 
Ti02 thin films. The diamagnetic contribution of the substrate, as well as that of 
the core electrons have been subtracted from both cycles. 



BE (inset of Fig. 4.12) that is formally assigned to a Ti^"^ oxidation state 



(3d^ configuration). This shoulder is absent in the stoichiometric Ti02 single 
crystal surface. In the present case, it is rather important to observe that, in 
spite of the N-doping, the Ti^"*" component, has virtually the same weight in 
both samples, and therefore, we can assume that the content of 3d^ magnetic 
ions is the same in both samples. 

In the N-doped sample, two components (N^ and N^) are present in 



the N Is core level region (Fig. 4.12, right panel). The peak at about 
BE=400 eV can be assigned either to molecular nitrogen(N2) bonded to 
surface defects or to N bonded to surface O sites''^. The intensity of 
component increases by collecting the spectra at grazing angles, i.e., in the 
more surface sensitive conditions for the XPS probe, and therefore this peak 
is ascribed to a surface contamination. In turn, the peak is ascribed to 
bulk N atoms incorporated during the thin film growth. An estimation of 
the N content was carried out by considering the area of the Ti 2p, N Is, 
and O Is photoemission peaks, suitably weighted by the photoemission cross 
section of each atomic level and the analyzer transmission. It is found that 
the N content in the N-doped film is about 
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Figure 4.12: Left panel: Ti 2p XPS spectrum from the undoped (a) and the 
N-doped (b) Ti02 thin film. These spectra are compared with that collected from 
a stoichiometric Ti02 (HO) single crystal surface (c). The inset shows the Ti 
2p3/2 core line region and the vertical arrow indicates electronic states that can be 
ascribed to Ti^^. Right panel: N Is photoemission spectrum at normal emission 
(gray dots) and at 50° off- normal emission (black dots). 
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The Ti L2,3-edge XAS spectrum is shown in Fig. 4.13, top panel. This 
spectrum is in agreement with those already reported for rutile with a set of 
bands that are related to the L3 edge (A, B, and C) and to the L2 edge (D 
and E). The A and D bands are ascribed to states arising from crystal 
field splitting, while the B, C, and E bands are ascribed to e^ states. The Pi 
and P2 prepeaks originate from multiplet splitting for the 2p^3d^ final state 
configuration. With respect to the undoped sample, the N-doped sample 
shows the same energy and intensity for all peaks but a larger width. The 
same remarks hold also for the O K-edge XAS spectra (Fig. 4.13, middle 
panel). Both spectra are well reproduced by DFT calculation of the undoped 
cell. The two spectra are identical but the peaks of the N-doped sample are 
slightly larger than those of the undoped sample, suggesting an higher degree 
of structural disorder. The fact that this extra broadening isn't observed in 
XPS spectra can be related to the different probing depth of XAS and XPS. 



Finally the N K-edge XAS spectrum is shown (Fig. 4.13, bottom panel). 
This spectrum is rather similar to that measured for oxygen, presenting a 
sequence of five bands (A to E). However a closer inspection shows that two 
additional features (Ai and Bi) can be detected and the C, D, and E peaks 
appear to be shifted with respect to the O Is XAS spectrum. DFT supercell 
calculations are able to consistently reproduce most of these features. In 
particular, the introduction of substitutional N determines a transfer of the 
spectral weight toward the low photon energies in the calculated A band and 
toward the high photon energies in the calculated B band, consistently with 
the observed Ai and Bi experimental features. Therefore, the Ai states can 
be regarded as N-related, donor, levels. For a very low concentration of Vqs 
the spectroscopic signature of the electronic states can be rather elusive and 
therefore one has to consider also ResPES techniques to probe their existence 
and determine the BE. Indeed, for both samples, the study of the occupied 
electronic states in the valence band has been carried out with ResPES with 
photon energies across the Ti L2,3-edge. 

In Figs. 4.14 a) and 4.14[ b), for each sample, we show the valence band 
photoemission spectra collected with photons below the Ti L2,3-edge (/iz/=450 
eV, off-resonance) and with photons corresponding to the A peak [hu =458.2 
eV, on-resonance) of the Ti XAS spectrum of (Fig. 4.13, top panel). As can 
be observed, in the resonance conditions three photoemission peaks (A, B, 
and C) results to be enhanced. Peak C is ascribed to the Auger emission, 
peak B to the resonant enhancement of Ti electronic states hybridized with 
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Figure 4.13: Top panel: Ti L2,3-edge XAS spectra collected from the undoped 
(a) and N-doped (b) rutile thin films. Middle panel: O K-edge XAS spectra col- 
lected from the undoped (a) and N-doped (b) rutile thin films, compared with the 
undoped cell DFT calculation. Bottom panel: N K-edge XAS spectrum collected 
from the N-doped rutile thin film, compared with the supercell DFT calculations. 
The resolution of the monochromator was set to 0.24 eV at the Ti L2,3-edge photon 
energy. 
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oxygen, and peak A to the contribution of Ti^^ induced by Vq. Tlie evidence 
of these in-gap states further supports the assignment of the tail on the 
low-BE side of the Ti 2p3/2 photoemission peak to Ti'^+ contribution. The 
N-doped sample shows a similar peak enhancement. 

The A features of both samples have been extracted by subtracting from 
the ResPES spectrum the off-resonance spectrum. The results are shown in 
Fig. 4.14[ c) (Vo-states). The two peaks present identical peak energies and 



intensities, providing an additional proof that the vacancy-induced states 
contribute to both samples with equal weight. The comparison of the off- 
resonance valence band spectra (Fig. 4.14[ c)) also allows to highlight small 



differences in the spectral weight that are well reproduced by DFT calcula- 
tions (Figs. 4.14K d) and 4.14[ e)). In particular, the doped sample displays a 



higher DOS with respect to the undoped sample for 2<BE<4 eV, and below 
10 eV. On the basis of the analysis of the DOS projected over the N-derived 
orbitals (Fig. 4.14K f)), the contribution at the top of the valence band, that 



adds electronic states (labeled as N-states) at low BE as compared to the 
pure rutile, is ascribed mainly to N. 

In summary, we have shown that N-doping of Ti02 rutile yields (i) an 
enhancement of ferromagnetic behavior, (ii) a replacement of about 4% O 
atoms in the rutile lattice with (substitutional) N atoms, and (iii) acceptor 
and donor levels, that determine a reduction in the rutile energy gap. This 
finding may have consequences on the enhancement of FM properties. In 
fact, the band gap reduction may favor the overlap of Vq states with the 
empty conduction band, rendering the latter magnetic. Moreover, a highly 
locahzed deep-gap defect orbital could lead to only a very-short-range defect- 
defect magnetic interaction that would fail to percolate through the sample. 
In turn, because of the band gap reduction, the Vo states can loose part of 
their deep-gap character and display longer-range interactions, creating more 
favorable conditions for the onset of magnetism. 

The different broadening of XAS spectra, as compared to the XPS one, 
suggest also that the defects could be localized away from the surface (or at 
least, away from the probing depth of photoemission technique). Moreover, 
as already shown in Chapter 2, the Vq states are deeply localized inside the 
bad-gap, while some theories for the d*^ ferromagnetism requires the creation 
of a impurity metallic band® that should be detected at the Fermi level. 
More details about the theoretical models are given in the next section. 
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Figure 4.14: Photoemission spectra collected from the undoped (a) and N-doped 
(b) samples in the valence band region. Spectra drawn with a thin line have been 
collected with photons below the Ti 2p-3d absorption threshold {hi' =455 eV), 
whereas spectra drawn with a thick line have been collected in resonance conditions 
(hi^ =458.2 eV). (c) Enlarged view of the off-resonance valence band spectra of the 
undoped (thin line) and N-doped (thick line) Ti02 rutile. (d) Calculated DOS for 
the undoped (thin line) and doped samples (thick line) samples, (e) Convolution 
of the calculated DOS with a Gaussian curve (FWHM=0.6 eV). (f) N-projected 
calculated DOS for N-Ti02. The calculated DOS curves have been shifted by 
about -3.5 eV to match the measured DOS. 
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4.6 Theoretical models 

The first attempt to explain the magnetism in DM0 were based on qual- 
itative considerations about the oxides electronic structure. Coey e al. ® 
originally proposed a model in which the magnetic impurities interact to- 
gether with the help of an impurity band created by Vos. This model is 
based on the assumption that the extra charge given by the oxygen vacancies 
is localized at the defect site and can be represented with an hydrogen-like 
atom. The radius of the Vo "wave-functions" then should be greatly en- 
hanced by the effect of the local dielectric constant, which is usually high in 
oxides. According to the authors, the FM should be present also when elec- 
trons remain localized in a narrow band because of correlation and potential 
fluctuations, thus justifying the magnetic interaction in insulating material 
(like most DMO's). 

Recently, the same authors® extended the model also for the (undoped) 
d° system, starting from considerations on the Fe doped Ti02 system. The 
magnetism should be of Stoner-type and takes place in the spin-split impurity 
band of percolating defects structure, with an half-metallic behavior. In this 
picture, only a 1% fraction of the volume should be involved in the FM. 
The Curie temperature value should be related to the splitting energy of 
the impurity band (in order of few mEv) and should be consistent with 
an high Tc ferromagnetism. Such kind of description fits well with our 
experimental findings, although is rather qualitative and not based on strict 
ab-initio calculations. 

A more systematical approach has been developed by J. Osorio-Guillen 
et al. for the undoped Hf02^ and CaOEB cases. The authors followed 
a stepped procedure: the first step is the individuation of the defects that 
can carry a magnetic moment; the second is the evaluation of the possible 
(both AF and FM) magnetic coupling through DFT simulations; the third 
is the evaluation of the percolation threshold for the desired lattice; the last 
is the evaluation of the lattice stability with the resulting density of defects. 
The authors identify the cationic vacancies as the only source of magnetic 
moments and estimate that the minimum density of defects needed for the 
percolation is around 5% for CaO and 13.5% for Hf02, thus much higher than 
what is allowed in an equilibrium-growth process. The latter conclusion is 
of course due to the choice of the cationic vacancies as the only magnetic 
moment carriers, while the most common defect in d° systems is the oxygen 
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vacancy. 

Summarizing, the theoretical models supports a defect-based origin of the 
FM. Because of the weak experimental saturation magnetization of typical 
d" systems, the oxygen (or cationic) vacancies should be spatially confined 
into a small region of the film volume. 



4.6.1 DFT on Ti02 supercells 

To better investigate the origin of magnetism in the rutile N-doped Ti02, 
another set of DFT calculations has been carried out in the spin-resolved 



LSDA-I-U® approximation, in the same supercell (see Fig. 4.15) described 
in Chapter 2. The purpose of these calculations is to further explore the 
contribution of nitrogen-doping to magnetism, apart from the band-gap re- 
duction. Four different supercell geometries have been considered: 

• an undoped cell; 

• a 2x2x2 cell with a substitutional N atom (labeled Nsub)', 

• a 2x2x2 cell with an oxygen vacancy (Vo); 

• a 2x2x2 cell with a substitutional N atom and an oxygen vacancy 

The same computational approach (PAW atomic description, with ABINIT''^ 
code) has been adopted; the cell geometry has been relaxed to the equilib- 
rium position in each case. A strong distortion of the lattice positions is 
expected, especially in the Vq case. The results for the electronic charge 
densities are shown in Fig. |4.16 with reference vectors that point out the 



relaxation effects induced by Vo and substitutional N. The triangle depicted 
in Fig. 4.16[ a) marks the three titanium atoms that surround an oxygen one 



in the rutile structure; the removal of this oxygen atom or its substitution 
with a nitrogen one directly influences the charge at the vertex of the triangle 
and ultimately the electronic structure. 

The relaxed geometric parameters of the different calculations are listed 
in Table 4.4 and are consistent with other results in literature'^. The Nsub 



supercell (Fig. 4.16[ b)) displays only a minor distortion, while the introduc- 



tion of the oxygen vacancy (Fig. 4.16 c)) clearly elongates the shape of the 
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Table 4.4: Structural parameters after the lattice relaxation. These parameters 



are labeled in Fig 4.16 the lengths are in A units. 



triangle, pushing the central atom away from the Vq. The maximum distor- 
tion is found in the case of the Vq plus the substitutional N: the central Ti 
atom is clearly shifted towards the nitrogen atom, which can accommodate 
the extra charge given by the Vq. This structural modification can be seen, 
for example, by looking at the a and d' values of the N^ub+Vo case in Table 



4.4, as compared to the other supercells. 




Figure 4.15: Supercell (2x2x2) geometry used for LSDA+U calculations. The 
single unit cell is shown in the lower left part. The supercell is divided in half by 
the (110) plane (in red) which contains the charge density slice shown in Fig. |4.16 



The total ionic charges of the Ti and N atoms near to the Vo have been 
roughly estimated with the Hirshfield method, as implemented in the CUT3D 
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tool included in the ABINIT package. As expected, in the undoped cell the O 
and Ti total charges are in a -1:2 ratio (0^~ , Ti'*+). In the Ngub supercell, the 
nitrogen atom keeps the same valence of the substituted oxygen (N^~). In the 
Vo supercell, the two electrons left from the Vq are distributed mostly in the 
surrounding three Ti atoms; the central Ti switches to an approximate Ti^+ 
ionic state, while the other two to a Ti^-^"*" state. Finally, in the Nsufe+Vo 
case the central Ti atom reverts to the Ti^"*" state, while its Vo electronic 
charge in excess is taken instead by the nitrogen, whose ionic state, evaluated 
by the Hirshfield method, is (i.e. with a filled 2p shell). 

The consequence of this charge redistribution on magnetism can be seen 
by looking at the spin-resolved DOS, shown in Fig. 4.17[ The spin-up and 
spin-down DOS of Vo supercell (see Fig. 4.17[ b)), like the undoped and the 
Ngufc, do not show any difference, reflecting an antiferromagnetic alignment 
of the spin on Ti atoms. The total magnetic moment of these configura- 
tions is thus zero; this also explains why non-clustered Vos have not been 
taken into consideration in theoretical calculations for the Hf02 and CaO 
cases''^ The two Vq"^ electrons magnetic moments are antialigned due 
to the exchange interaction. In fact, as suggested also by Coey and collab- 
orators'^', only clustered Vos can display magnetism. A non-zero magnetic 
moment is expected in neutral (Vq) or single-charged (V^") vacancies''^. 

However, a nitrogen atom near the Vo can attract one of the Vq"^ elec- 
trons, leading to a net magnetic moment. The spin-resolved DOS on the 
NsMfe+Vo cell (Fig. 4.17K c)) displays both a band-gap reduction (0.2 eV) and 
spin-polarized in-gap states. The spin-difference charge density plot shown 
in Fig. 4.18 shows how the presence of nitrogen can quench the spin-down 
density (blue color in Fig. 4.18[ a)) on the central Ti atom, thus leading to a 
non-zero total magnetic moment. 

The nitrogen doping on rutile Ti02 can increase the FM not only through 
the reduction of the band-gap, that may favor a long range interaction be- 
tween Vos. In fact, the substitutional N atoms can host part of the Vo 
excess charge, increasing the number of magnetically active sites. Since ev- 
ery oxygen atom is bonded to three titanium atoms (and thus is shared by 
three oxygen octahedra), in a rutile lattice with a 4% of uniform substitu- 
tional N-doping the probability for a Vo to be near a dopant N is about 
60%. Therefore, in our films at least half of the Vos, even if isolated, should 
display an enhanced magnetic moment. In a film with a 6.6 % of uniform 
N-doping, the presence of a nitrogen atoms in the proximity of each Vqs is 
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Figure 4.16: Total electronic charge density along the (110) plane (see Fig. 



4.15) for the various LSDA+U calculations. The data on the undoped rutile(a) 
cell have been extended to match the 2x2x2 range for comparison. The values for 



the inter-atomic lengths labeled with the reference vectors are given in Table 4.4 
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Figure 4.17: Spin resolved DOS for the Vo (b) ^sub+^o (c), compared to the 
single rutile cell (a). The zero of the energy scale is set to the top of the valence 
band, without considering the in-gap defect levels. 

assured. This kind of mechanism is rather general and could be applied also 
to other N-doped magnetic oxides. 

4.7 Conclusions 

In conclusion, in this Chapter a review of magnetism in Ti02-<5 and 
N:Ti02-5 grown with RF-sputtering is given. A weak room-temperature 
ferromagnetism has been detected in each sample, despite the absence of 
magnetic ions (TM or RE) doping. In-vacuum annealing systematically in- 
creases the saturation magnetization, while the annealing in oxygen results 
in the opposite effect. This can be regarded as a proof of the Vo fundamental 
role in FM, also excluding the possibility of a contamination from the sample 
holder and the annealing stage. 

The weak intensity of FM cannot be explained with a uniform distribution 
of magnetic spins; FM should be then explained by a clustering of defects 
in particular areas of the samples, as suggested from theoretical works in 
literature. The saturation magnetization per unit volume is higher in thin 
films than in the corresponding thick ones, suggesting the presence of plane 
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Figure 4.18: Spin-difference electronic charge density along the (110) plane (see 
Fig. 4.15) for the Yq (a) N^^^+Vo (b) configurations. Blue color corresponds to 
negative values, red to positive one. 
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of clustered defects, either close to the surface or to the interface. XPS 
measurements on Ti 2p core levels, a surface sensitive technique, did not 
unambiguously display Ti^+ (or Vo) signatures, thus supporting the hypoth- 
esis of a buried FM interface defect layer in our films. RcsPES at Ti L-edge, 
taken on a rutile thin film, has proved to be an effective tool to underline 
the presence of electronic 3d^ states in the gap, even for a very small concen- 
tration of Vos; however, the observed Ti^+ in-gap states lie away from the 
Fermi energy and thus can be ascribed to isolated oxygen vacancies with a 
strongly localized electronic structure. DFT-I-U calculations on rutile super- 
cells further support this description. 

Upon nitrogen doping, the ferromagnetism in rutile Ti02-5 is greatly en- 
hanced. The main effect of the N-doping in the electronic structure is the 
band-gap reduction, which can be easily observed with optical transmission 
experiments. XPS and ResPES have been performed in order to underline the 
nitrogen related electronic structure in core-levels and in the valence band. 
The structural disorder induced by the N-doping can be seen through the 
broadening of XAS and XPS spectra. Finally, spin-resolved LDA-I-U calcu- 
lations suggest that substitutional N-doping can induce a magnetic moment 
in otherwise non-magnetic isolated Vq" through Ti oxidation. 

In conclusion, even if the magnetism in undoped Ti02-5 is still far from 
being exploited in future applications, it represents a challenging issue for 
solid state physicists. In order to understand the magnetism in Ti02-5, 
further experimental investigations are required, in particular to map the 
bulk Vo distribution. From the theoretical side, ab-initio calculations on 
large defects clusters are needed in order to give a realistic description of 
these systems. 
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Chapter 5 



The LaAlOs-SrTiOs 
heterostructure 



5.1 Introduction 

Lanthanum aluminate (LaAlOa, LAO in short) and strontium titanate 
(SrTiOs, STO) are ionic insulators, with a band gap respectively of 5.4 and 

3.2 eV. These materials belong to the perovskite group, which share the 
same chemical formula (ABO3) and a similar cubic crystal structure (see 



Fig. 5.1). These compounds are very interesting due to their relative simple 
crystal structure and to the possibility to have rare earths or transition metals 
at the cationic A and B sites, displaying many different functional proper- 
ties (metal-to insulator transitions, ferroelectricity, magneto-resistance and 
so on). Moreover, it is possible to use a perovskite as an epitaxial template 
for the growth of another one, creating a seamless amount of possible het- 
erostructures. LAO and STO are formally band insulators, since they are 
closed-shell compounds (4f^ for LAO and 3d° for STO) and are easily grown 
with chemical vapor deposition (CVD) or pulsed laser deposition (PLD); in 
fact, high quality STO samples are commonly used as a starting substrate 



for epitaxial growth. As shown in Fig. |5.1[ the LAO-STO interface can be 
p-type when the bulk STO is terminated with a SrO plane (hole doping) and 
n-type with a Ti02 plane (electron doping). In the latter case, an interesting 
phenomenon occurs: the interface of LAO-STO becomes conducting''^ and 
forms a quasi-2D electron gas (2DEG). The proof of the 2D character of the 
conductive layer was given later by means of the anomalous Shubnikov de 
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Figure 5.1: A schematic illustration of LaAlOa and SrTiOa crystalline and elec- 
tronic structure. 



Haas effect G3. 

Moreover, there's a thickness dependence of the metal transition: the 
2DEG is observed only when the capping is at least 4 unit cell thick '2^. The 
main difference between LAO and STO resides in the layer charge polarity''^: 
looking at the (001) planes (see Fig. 5.1), STO is a non-polar solid, since 
both Sr^"*" 0^~ and Ti^"*" Ol~ are charge- neutral, while LAO is a polar solid, 
as it is composed of La^+ 0^~ and Al^"*" Ol~ charged layers. The p-type 
interface is thus formed by SrO - (A102)^~ planes, n-type by Ti02 - (LaO)^"*" 
ones. The observed conductivity should be the response of the system to the 
diverging potential (the so called "polar catastrophe") created by LAO, both 
in the electronic and structural properties. 

In non-oxide semiconductors, the relaxation of polar discontinuity in het- 
erointerfaces is usually achieved by the so called "atomic reconstruction" 
process''^: the interface stoichiometry is altered by interdiffusion, point de- 
fects, dislocations and in general by a structural roughening. In oxides, the 
possibility of multiple valence ions allow also an "electronic reconstruction" 
that, in this case, should move electrons from the surface to the empty Ti 
d-level, creating 3d^ electronic states. A 2D lattice of electrons in a corre- 
lated material can raise phenomena''^ like MIT transitions, localized mag- 
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netic moments and even superconductivity; most of tliis phases have been 
observed, but not simultaneously in the same interface. 

In principle, both atomic and electronic reconstruction could be present 
in the LAO-STO case. There are experimental hints that support either 
one theory or the other. For example, there are many experimental proofs 
of interdiff usion ''^ (lanthanum ions moving inside STO, a form of atomic 
reconstruction), but this alone cannot be the solely responsible of the con- 
ductivity, since in principle should be present also in the n-type interface. 
Moreover, standard "polar catastrophe" (an electronic reconstruction pic- 
ture) alone can't explain the lack of conductivity for thickness below 4 unit 
cell threshold. The thickness dependence could be explained by a band bend- 
ing effect, induced by polarity discontinuity: the LAO VB DOS should be 
shifted to higher binding energies till, for a capping equal or major of 4 
u.c, the VB maximum is superimposed to the buried empty levels of STO. 
The conduction should be caused by a tunneling effect from surface to the 
interface. However, a band bending has never been observed in therm of 
core-level shift, while a shift of 3.2 eV (needed to spam the electronic gap in 
STO) should be easily observed. 

Finally, the sample quality deeply affects the transport measurement; 
an oxygen-poor growth atmosphere can induce oxygen vacancies and thus 
a 3D conductivity, while an excessively-rich one can even result in a 3D 
growth and thus in a different kind of heterostructure. In Nb-doped STO, the 
presence of oxygen-vacancies at the interface can induce superconductivity 
with Tc below 400 mK, in a similar way to the one measured on LAO-STO''^. 
Three phases are usually identified, according to the O2 partial pressure 
during growth (Fig. |5.2 ): one dominated by oxygen vacancies contribution 
(P02 — 10^^ mBar), one displaying superconductivity {P02 — 10~^ mBar) 
and one displaying a magnetic behavior (P02 — 10^'^ mBar). It is quite a 
challenging task to find a unified description of all this phenomena. 

Since the conduction is thought to occur at the interface, the filling of 
electronic states at Fermi level can't be observed with standard photoemis- 
sion techniques: in conductive samples, even in a synchrotron facility. X-ray 
photoemission might not show the Fermi level, being more sensitive to the 
uppermost insulating oxide. The problem can be solved, as in the Ti02 case, 
using the enhancement of selected spectral weight at resonance condition. 

In this Chapter, an analysis of photoemission data of insulating and con- 
ductive LAO-STO film is given; samples have been grown with PLD in a 
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Figure 5.2: Temperature depen- 
dence of the sheet resistance for 
SrTiOa/LaAlOs conducting inter- 
faces, grown at various partial oxy- 
gen pressures. Three regimes can 
be distinguished: low pressures lead 
to oxygen vacancies, samples grown 
at high pressures show magnetism, 
whereas samples grown in the inter- 
mediate regime show superconduc- 
tivity. Adapted from Ref. ^ . 

Po2=l"10~^ mBar oxygen partial pressure at the University of Twente by 
A. Brinkmann research group. In the first section an XPS analysis of the 
stoichiometry of LAO-STO as compared to single-crystal LAO and STO is 
shown; the second section is devoted to ResPES, XAS and XLD analysis; 
the third section is devoted to intermixing analysis by means of AR-XPS. 
Finally, the last section is dedicated to conclusions and future improvement 
in research. 

5.2 XPS data 

The XPS measurements have been carried out on the 3 u.c. and 5 u.c. 
LAO-STO samples, plus two reference single- crystal LAO and STO samples 
terminated with the (001) surface. In addition, an insulating 5 u.c. sample, 
grown at 10^^ mBar O2 partial pressure, has been analyzed. XPS has been 
used to measure the core-level electronic structure and to evaluate film sto- 
ichiometry. The STO and LAO samples have been treated with Argon ion 
sputtering in order to see the effect of oxygen loss. The XPS data shown 
in this section have been collected at Surface Science Lab of the Univer- 
sita Cattolica Brescia, with a non-monochromatized dual-anode PsP x-ray 
source; the Mg ka line (/iz/=1256.6 eV) has been used to achieve a better 
resolution (about 0.7 eV), while the Al line (/iz/=1486.6 eV, resolution of 
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Figure 5.3: XPS survey spectra taken in single-crystal LAO, STO and on the 3 
u.c and 5 u.c. LAO-STO samples. Peaks are labeled on the figure. Spectra have 
been collected with a Mg k^ x-ray source and have been normalized to the oxygen 
Is peak. 



0.8 eV) has been used when the maximum probing depth was needed (mainly 
for AR-XPS). Charging effects have not been detected during XPS, since this 
kind of X-ray source are rather unfocused. The analyzer was a SCIENTA 
R3000 XPS/UPS/ARPES, with a maximum energy resolution of 25 meV. 
For XPS, this analyzer has been used in the "transmission mode", which 
maximize the transmittance and work with a 30° acceptance angle; "angular 
modes" are also available, and have been used for ARXPS. A survey view of 
the detectable core levels in 3 u.c. and 5 u.c. samples is shown in Fig. 5.3[ 

Carbon Is peak is due to the usual atmospheric contamination: samples 
hav not been cleaned up since vacuum annealing or Ar+ sputtering might 
have changed the oxygen stoichiometry. Moreover, C Is peak can be used as 
reference for the binding energy scale (C Is binding energy for carbon con- 
tamination is fixed at 284.8 eV). Even in the survey spectra. La 3d and 4d 
peaks show already a complex structure which is not due to different oxida- 
tion state but to well-known multiplet structure ''^S' . The different intensity 
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Figure 5.4: XPS survey spectra taken on Sr 3d, Al 2s and La 4d on LAO and 
heterostructures. Spectra have been normalized to the AI 2s area to allow a better 
comparison of lanthanum over aluminum ratio, which are given in the legend. 
Calculated ratios are already normalized with the photoionization cross section. 
Experimental condition (i.e. the total resolution) has been kept the same for all 
the samples. 



of LAO vs STO structures in the two samples is already visible in Fig. 5.3 



basically, in the 5 u.c. sample Ti and Sr peaks are more attenuated with 
respect to the Al and La ones due the larger LAO capping. 

Since the 80-140 eV (in BE) region contains, in a relatively short range, 
most of the core level peaks required for a quantitative analysis, high resolu- 
tion spectra (see Fig. 5.4) have been collected on this area to carry out the 
element quantification. The peak areas have been measured after a Shirley- 
type background subtraction and have been normalized according to the 
photoionization cross sections''®. First of all. Fig. 5^ spectra clearly show 
that the insulating 5 u.c. is very different from the other samples because, as 
already pointed out in the introduction, its growth regime is 3D instead of 
layer-by-layer. This sample still retains a crystalline order, as will be shown 
with AR-XPS data in the next section. The Al/La ratio in heterostructures is 
different from the LAO crystal one, showing a lanthanum deficiency that can 
suggest a La interdiffusion; anyway, LAO-STO samples are terminated with 
AIO2 plane while single-crystal LAO is usually terminated with a mixture of 
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Figure 5.5: High-resolution Ti 2p XPS taken on the 3 u.c LAO-STO, conductive 
5 u.c. LAO-STO and untreated single-crystal STO. The black arrows identify the 
position of Ti'^"*" states, increasing from STO to the 5 u.c. sample. Similar data 
have been published in Ref . *^ . 



AIO2 and LaO phases, resulting an higher (closer to 1:1) ratio. The measured 
value for LAO-STO is compatible with the one that can be calculated with 
the AR-XPS model (see next section) and thus is not a sufficient proof of 
selective lanthanum diffusion. Moreover, data in Fig. 5.4 have been collected 
in the "transmission" mode, thus the results are the average of the possible 
diffractive effect in the analyzer acceptance angle. In order to have pre- 
cise theoretical results, one should rely on precise photoelectron-diffraction 
calculation to verify peak intensities enhancement at normal emission. 

Probing the Ti/Sr ratio can be a more difficult task for two reasons: first, 
the most prominent peak (Sr 3d and Ti 2p) are 300 eV apart and the analyzer 
transmission function could not be constant; second, the total area of Ti 
2p peaks has to be measured carefully accounting for the charge-transfer 
satellites. The results for this two peaks ratio are given in the AR-XPS 
section. 



In Fig. |5.5| the Ti 2p XPS peaks of the 3 u.c, 5u.c. and pure STO samples 
are shown. The titanium 2p peaks are almost identical to that expected for 
a Ti^"*" ion. However, a weak spectral feature (indicated by the black arrow) 
is detectable on the lowest BE side of 5 u.c. spectrum in the position usually 
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associated to Ti^"*" electronic states. These peaks can be detected only by 
comparison with different samples and can be easily confused as an addi- 
tional experimental broadening of Ti 2p3/2 peak. Recently, these electronic 
states have also been measured with hard X-rays (HAXPES) experiment to 
evaluate the Ti interface reconstruction in LaAlOa-SrTiOs ''^ with a large 
LAO capping. 

In order to underline the effect of an oxygen loss in XPS spectra, in Fig. 



5.6 data taken on the clean STO after subsequent Ar"*" IkV sputtering treat- 



ment are shown. The main effects on peak shapes are indicated by the black 
arrows: Ti 2p3/2 gains the low BE component, Sr peaks become broader and 
O Is peaks lose the component at higher BE due to surface contamination. 
The spectra (collected at the same resolution) of heterostructures are also 
shown; in particular, the Sr peaks are not as sharp as pristine STO and are 
comparable to the spectrum taken after one sputtering session. This can 
indicate a sort of disorder (cationic exchange or oxygen vacancies) around 
strontium atoms at the interface; other indications supporting this picture 
are given in the ARXPS section. 

A similar analysis has been carried out on a clean LAO (001) sample; 
results are shown in Fig. 5.7, As expected the intensity of C Is peak, which 
is associated to adventitious contamination, reduces upon sputtering; the 
shoulder of O Is is quenched as in the STO, but the overall width of main 
peak also increases due to the increased disorder. The La 4d peak width (in a 
similar way of Sr 3d in STO) gets broader but surprisingly the La 36.^/2 seems 
to behave in the opposite way. In fact, the splitting of La 3d in Fig. 5.7 is not 
due to spin-orbit interaction but it is caused by multiplet structure, common 
in lanthanum (and in the other RE) oxides. The different broadening in this 
case can be explained by the presence of metallic lanthanum, whose La 3d5/2 
peak-shape is given only by a single peak that raise the low BE side of the 
overall XPS spectrum. Anyway, the shape difference in LaAlOa core level 
can give little information to the LAO-STO case. 

Finally in Fig. 5.8, 5.9, 5.9 the valence band (VB) spectra, taken with 
the Mg ka X-ray source, are shown. The shallow core levels are labeled in 
Fig. 5^; the O 2s spectral weight is located mostly at Si peak, Sr 4p is an 
unresolved doublet below S2 peak and La 5p is splitted in S2 and S3 peaks. It 
is possible to describe the LAO-STO spectra as the linear combination of the 
single-crystal data; the spectra of LaAlOs and SrTiOs have been normalized 
in intensity, shifted and summed in order to minimize the difference with the 
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Figure 5.6: Ti 2p3/2 (a), Sr 3d (b) and O Is (c) XPS spectra measured after 
subsequent Ar"*" sputtering treatments on the STO sample, as compared to LAO- 
STO data. The black arrows mark the changes due to sputtering effects. LAO- 
STO samples have not been sputtered in order to avoid structural damages or the 
creation of oxygen vacancies. 
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Figure 5.7: XPS spectra from selected core level spectra taken during subsequent 
sputtering session on a single crystal LaAlOs sample. The colors of experimental 
data range from black (not treated) to grey (heavily sputtered). The black arrows 
underline the modification caused by sputtering. 
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Figure 5.8: XPS valence band spectra on single crystal (c and d) and heterostruc- 
tures (a and b). Both the VB and the shallow core levels are shown. 
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multilayers. The results of this "fitting" procedure are given in Fig. |5.9[ there 
is a minimal shift difference of 0.2 eV in the two linear combinations, which 
should reflect the difference in the band bending between the 3 u.c. and the 
5 u.c. samples. Such a small value, not enough to justify the conductivity, 
is similar to the one that can be found with ResPES measurement, and 
will be discussed in the conclusion. The LAO/STO ratios indicated in Fig. 



5.9 are useless by themselves, since the absolute XPS intensity can vary in 



different spectra due to the different thickness of the LAO capping; aniway. 



it is possible to evaluate the difference in thickness by using the formula 2.4 
The STO signal attenuation of 5 u.c. respect to the 3 u.c. at normal emission 
can be written as: 

-J = e ^ (5.1) 

J-STOZu.c. 

The fitting procedure in this case gives a ratio of about 0.59. For electrons 
with a kinetic energy of 1250 eV, the TPP-2M formula in pure LAO gives 



an inelastic mean free path A = 23.32 A. By inverting Eq. 5J^, a difference 
in thickness oi z = 12.30 A is found, much higher than the expected 2 u.c. 
(7.6-7.8 A). This indicate that signal of the 5 u.c. sample is much more 
attenuated than expected respect to the 3 u.c. case. 

In fact, the problem here is not related to the samples quality, but to 
the oversimplification of Eq. |5.1 which doesn't take into account the elastic 



scattering. By replacing A with the calculated EAL, as described in Chapter 
2, an extra-thickness of nearly 9 A is found, in good agreement with the 
expected value. 



In Fig. 5.10 the experimental LAO and STO valence band spectra, com- 
pared to the total density of state calculated with DFT on the cubic per- 
ovskite single cells, are shown. An LDA-I-U approach has been used in order 
to recover the correct band g implemented in the ABINIT® code. 

The parameters used in the calculation are® U5to=8.0 eV, J5to=0.8 eV, 
Ulao= 11-0 eV and Jlao=0.68 eV. Since the VB is mostly composed by just 
one component (oxygen 2p states), the cross section of different electronic 
levels has not been considered. The agreement with calculation is fairly 
good; LAO states in particular are marked by the extra spectral weight at 
BE=11 eV, which is found also in DFT calculation. This good agreement 
can help to achieve a better energy alinement of the separate LAO and STO 
contributions in LAO-STO, especially in the ResPES measurement. 
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Figure 5.9: 3 u.c. and 5 u.c. LAO-STO valence band fitted with a linear combi- 
nation of LAO and STO spectra. The starting intensity of the single components 
has been normalized to the total VB area and a maximum relative shift of ±1 eV 
has been considered in the fitting procedure. 




15 10 5 
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Figure 5.10: Single crystal LaAlOs and SrTiOs valence band spectra, compared 
to DFT calculations of the density of states. The calculations have been done on 
a relaxed cubic perovskite single cell with the LDA+U approach. 
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Figure 5.11: 3 u.c. and 5 u.c. LAO-STO Ti L2,3 edge absorption spectra, 
compared with multiplet calculation. XAS spectra have been recorded at 30° 



incidence with linear polarized light (horizontal and vertical, see Fig. 5.14), and 
calculations have been done for an isolated Ti^"*" in a slightly distorted octahedral 
(D4/1) symmetry, following the work of Salluzzo et. al.'*"'^'. X-ray linear dichroism 



results are given in Fig. 5.12 



5.3 XAS and XLD 

X-ray absorption spectra and X-ray linear dichroism spectra have been 
collected at the BACH beamline of the ELETTRA synchrotron. Measure- 



ments have been carried out in the same configuration of XPS (Fig. 5.14), 
with a photon resolution of 0.15 eV at hi'=455.0 eV, at room-temperature. 
XAS spectra have been taken with a completely in-plane s polarization (verti- 
cal polarization) and with a partially out-of-plane p polarization (horizontal 
polarization). Total yield (i.e. drain current) has been the most reliable 
XAS detection mode, given the insulating nature of 3 u.c. samples. In Fig. 



5.11 XAS results are shown on Ti L2,3 edge for 3 u.c. and 5 u.c. LAO-STO 



samples. 



As shown in Fig. 5.11, the shape of the Ti L2,3 absorption edge can be 



well reproduced by a multiplet calculation for a Ti^"*" ion in an 0/j symmetry. 
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with 10Dq=1.6 eV. The main difference from Ti02 spectra is the absence of 
L3 Eg peak sphtting (see Fig. 2.8), due to the higher symmetry of the [TiOe] 
cluster in the heterostructure. Salluzzo et ah ^ have already reported a 
very weak dichroic signal at LAO-STO Ti L-edge, as a difference of circu- 
larly and linearly polarized XAS. Dichroism is not expected in an perfect Oh 
environment, so that theoretical spectra (Fig. 5.12) have to be calculated by 
slightly lowering the symmetry group to a distorted octahedral (D4/i), with 
ds=dt=25 meV. However, from an experimental point of view, such differ- 
ences are so small that could be generated by a minor shift (less than 0.1 
eV) of photon energy during data acquisition. 

In fact, the difference of two identical spectra shifted by a small quantity 
on the X-axis could give an effective dichroic signal, caused by the experimen- 
tal set-up and not to by a real physics effect. In a genuine dichroic effect, the 
spectral shape has to be different (and not only translated) on in-plane and 
out-of-plane polarization. A check of the linear dichroism can thus be eval- 
uated by considering the difference between the s-polarized and the shifted 
p-polarized spectra, till a minimum of total intensity is found as a function 



of displacement. The results of this procedure are given in Fig. 5.12 



The subtraction of Fig. 5.11 XAS spectra yields a strong measurable XLD 
signal for both samples, but a careful check (the results of minimization are 
labeled as displaced spectra of Fig. 5.12) reveals that in the 5 u.c. the XLD 
signal persist even after the displacement, while in the 3 u.c. the minimum 
difference spectra is only positive, indicating the absence of a real shape 
differences in the p and s spectra. 

These findings could indicate that the extra-energy caused by the disconti- 
nuity is relaxed through a ferroelectric-like distortion "ED of the Ti octahedron 
at the interface, and that this distortion is higher in conductive sample than 
in the insulating one. Moreover, both XAS and XLD data don't show any 
trace of Ti^+ related features (theoretical spectra can be always calculated 
with a Ti^"*" configuration only). Since the XPS data provide evidence for 
Ti^+ states, the differences in results can be ascribed to the different probing 
depth of the two technique, the XPS being the most sensitive to interface. 
The octahedron distortion can be a consequence of the strain effect due to 
the lattice mismatch''^. 



In Fig. 5.13 the XAS spectrum at lanthanum M4^5 edge of 3 and 5 u.c. 
samples is given. These spectra are composed by of a simple 3 peak structure, 
two main peaks (3d3/2 and 3d5/2) and a weak small pre-peak. The shape can 
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Figure 5.12: 3 u.c. and 5 u.c. LAO-STO Ti L2,3 edge XLD spectra, compared 
to multiplet calculations. First, the spectra of Fig. 5.11 have been subtracted 



from each other as measured (dichroism spectra); then a relative shift has been 
imposed in order to minimize the XLD signal ('effective' dichroism spectra). The 
conductive 5 u.c. sample shows a detectable dichroism, while minimization on 3 
u.c. spectrum give a non dichroic signal (always positive) indicating that the two 
s and p original spectra are identical in shape. Multiplet calculation for a Ti^"*" 
ions in a symmetry are shown in the bottom. 
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be easily calculated with multiplet calculation''^, and doesn't correspond to 
the presence of metallic lanthanum. No dichroic signal has been detected. 
Oxygen k-edge XAS is not reported here since it contains contribution of 
STO, LAO and the surface contamination. Moreover, charging effect on 
insulating oxides are enhanced at oxygen absorption edge, preventing the 
collection of the 3 u.c. spectra. 



5.4 ResPES 

The ResPES data have been collected in two different sessions at BACH 
beam line with two sets of identical LAO-STO samples; results of first beam- 
time (already published in Ref. ''^) have been used to better focus on the 
data collection of the second beam time. Experimental results of the two 
beam-times are given in two separated paragraphs. To achieve the maximum 
"bulk" sensitivity during measurement, photoelectron have been collected at 



normal emission angle (see Fig. 5.14), thus forming a 30° angle with the soft 
X-ray beam. 



a: 3d view b: Top view 

z 




h. pol. 



detector 

Figure 5.14: Analyzer and sample geometry during photoemission and XLD 
measurement; in A the 3D view, in B the top view. X-ray vertical and horizontal 
polarization (v. pol. and h.pol.) are marked by green and blue arrows. 



This set-up, which is the same used for XLD measurements, is suitable 
for measuring dichroic effects in this sample; in particular, the "vertical" 
polarization (Fig. 5.14[ b)) has an electric field that oscillate entirely on 
sample plane, while the "horizontal" polarization has a strong out-of-plane 
component. One of the aim of this work has been to find possible dichroic 
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Figure 5.15: 3 u.c. (a) and 5 u.c. (b) LAO-STO resonant photoemission spectra 
taken at Ti L edge. Spectra contains shallow core levels (around BE=22 eV) and 
valence band (around BE=7 eV). Main (Ti^"*") resonances are marked with the 
letter B. Taken from Ref.C^. 



effects in photoemission, due to the quasi-2D nature of the conductive layer. 



5.4.1 First data set 

As in the Ti02 case, ResPES have been collected at the Ti L2 3 edge. The 



results for both the 5 u.c. and 3 u.c. LAO-STO are shown in Fig. 5.15 



The ResPES spectra Fig. |5.15| show a resonant behavior of the peak at 
BE=6 eV (marked as B in both ResPES set of data), which is ascribed to Ti 
states hybridizing with 0-states in the valence band, as described in Chapter 
3 for the Ti02 case. The resonance is most intense for the 3 u.c. sample (Fig. 
5.15[ b)), where the Ti is closer to the surface than for the 5 u.c. sample. The 



ResPES also reveals in-gap electronic states marked as A in Fig. 5.16 The 



analysis of these states is accomplished by examining two spectra collected 



with a photon energy of 457 eV, as shown in Fig. 5.16 At this excitation 
energy, only the Ti^"*" seems to resonate. 

Apparently, the in-gap states have the same weight in both samples. How- 
ever, because of the thicker LAO overlayer in the 5 u.c. sample, the effective 
density of the in-gap states is larger for the conducting 5 u.c. sample. The 
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Figure 5.16: Valence band spectra of the 3 and 5 u.c. samples collected at a 
photon energy of 457 eV. The inset shows the in-gap states for both samples. 
Taken from Ref.C3. 
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Figure 5.17: XAS of the 5 u.c. sample (thick Hne) and of LaTiOs (adapted 
form Ref. 'ED, triangles). CIS profile of the A and B features for the 3 u.c. sample 
(empty and filled diamonds, respectively). CIS profile of the A and B features for 
the 5 u.c. sample (empty and filled circles, respectively. Taken from Ref. '1^. 



excess of in-gap electronic states is consistent with the Ti'^+ states found in 
XPS Ti 2p core level data of Fig. 5.5 If we assume that the only difference 
between the two samples is the thickness of the LAO overlayer, the attenu- 
ation of the signal based on the estimated for electrons with 457 eV kinetic 
energy would give an j^^.c. ^ j5u.c. _ -^ gj ratio between the intensities I of the 
A peaks of the two samples. Thus, as the measured intensities appear to be 
the same, the density of in-gap states for the 5 u.c. sample results to be 1.97 
times that of the 3 u.c. sample. These estimates will be improved with the 



second data set results. Finally, CIS spectra are given in Fig. 5.17 



The CIS spectra are shown for the in-gap states (Fig. 5.17, A feature), 
as well as for the most resonating feature at BE=6 eV (Fig. 5.17, B feature). 
The peak at 458.2 eV is ascribed to t2g states arising from crystal field split- 
ting, while the broad peak at 460.6 eV is ascribed to eg states. As already 
explained in Chapter 2, the CIS spectra of the B feature follow the intensity 
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of the Ti La-edge XAS (only the XAS for the 5 u.c. sample is shown, the 
other being quite similar) in agreement with their origin from Ti^"*" electronic 
states in the valence band. As in the case of defective Ti02, the CIS spectrum 
of the A band shows a different behavior, with a quenching of the intensity 
below the t2g peak and a shift toward lower photon energies of the e^ band. 
In spite of the weaker signal with respect to the CIS spectrum extracted form 
peak B, the CIS profiles of peak A show a higher intensity where the most 
intense XAS bands of a Ti"*"^ ion are found (for comparison, the Ti L3 edge 
of LaTiOa is given. Fig. 5.17, thick line). This adds further evidence to the 
3d^ character of the in-gap states. However, the experimental noise on CIS 
spectra of Fig. 5.17 doesn't allow any prediction about the symmetry of Ti 
bonding at the interface. 



5.4.2 Second data set 

In order to improve the quality of CIS and in-gap state spectra another 
set of ResPES data has been collected; in this second beamtime also the 
polarization properties of VB photoemission have been proved. Most of the 
experimental data (CIS and ResPES spectra) shown in this paragraph have 
been taken with in-plane polarized X-rays, which give the maximum intensity 
on Ti resonances. Further comment on linear polarization effect will be given 
at the end of the paragraph. Following the case of Ti02, the Ti L-edge 
ResPES more spectra have been collected in the pre-edge region where only 
the Ti'^"'" states should be present. The X-ray flux intensity has been reduced 
for insulating samples in order to control the charging effects, resulting in an 
higher noise on experimental data. ResPES data for the 5 u.c. sample are 



shown in Fig. 5.18 



The maxima of in-gap state intensity, marked with black arrow in Fig. 



5.18 (a), are at different photon energy than main resonance peak. In order 



to better visualize the in-gap states, in Fig. 5.18[ (b) an enlargement of 



ResPES image is depicted, with a different color scale. Panel (c) of Fig. 



5.18 depicts selected photoemission spectra from panel (a) data, measured 
in out-of-resonance conditions (/;,z/=452 eV) and at the maximum intensity 
of Ti^"*" resonance (/iz/=458.1 eV) and Ti^"^ (/ii/=459.6 eV); the horizontal 
dashed lines mark the single spectra position on the main image. Peaks 
Ai,A2 and A3 are the shallow core levels, B and D are the most prominent 
VB features, C is the position of Ti^"*" resonance peak maximum and E is the 
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Figure 5.18: Second ResPES data-set at L-edge on 5 u.c. LAO-STO sample. In 
panel (a) ResPES data in a 2D picture and in panel (b) a detailed view on 

Ti3+ 

defects; on the left axis XAS (black line) is depicted. In panel (c), spectra taken 
at out-of resonance (blue line, blue dashed horizontal line of panel (a)), maximum 
of Ti^^ spectra (red line) and maximum of Ti"^^ (black line) conditions. In panel 
(d) the detail of in-gap state spectra. 



Ti^"*" electronic state. Remarkably, the in-gap state appear to be at Fermi 
level, as expected from a conductive sample. 



As shown in Fig. 5.19 the ResPES of the 3 u.c. sample displays the 
same identical features of the 5 u.c, except for the relative intensity of main 
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VB resonance and shallow core levels; the peak marked by F in Fig. 5.19 



can be ascribed, as in the case of Ti02, to the LVV Auger emission and is 
more pronounced due to the lower LAO capping. Again, the defect peak 
is detected near the Fermi level, its intensity being the only difference with 
respect to the 5 u.c. case. 



In Fig. 5.20 the CIS spectra for 3 u.c and 5 u.c. are shown and compared 
to Ti L-edge XAS spectra on samples and on a reference LaTiOs crystal. The 



Ti^+ (Ti3+) CIS has been collected at C (E) peak of Fig. 5.18 (c) and Fig. 



5.19,(c). The similarities between LaTiOs XAS and Ti^+ CIS are evident, 
suggesting that of La-Sr intermixing may occur at the interface, yielding Ti^+ 
states as in Sr doped LaTiOa. Since the in-gap state intensity is lower in the 
3 u.c. sample, different scaling factors have been applied to CIS spectra of 
Ti3+. 

A first comparison between the 3 u.c. and 5 u.c. LAO-STO samples is 
given in Fig. 5.21 black line spectra have been taken in out-of resonance 
conditions {hu = 450 eV), blue line have been taken at the maximum of Ti^+ 
resonance {hu = 459.6 eV) and red spectra show only the resonating part 
given by the difference of the previous two. A normalization on the shallow 
core level areas has been done before the subtraction. By looking at the 
non-subtracted resonating spectra (black lines) of Fig. 5.21 it seems that the 
intensity of defect states is similar in both samples; in fact the Ti'^+ states 
have an effective intensity which is much lower in the 3 u.c. LAO-STO, since 
the titanium photoemission signal should be more attenuated in the 5 u.c. 
sample because of the larger capping. The calculated ratio is 
nearly 2:1, considering the IMFP only. 

This value can be affected by the relative normalization of the 5 u.c. and 
3 u.c. ResPES data; a better comparison can be obtained by looking at the 
resonating part only (green spectra of Fig. 5.21) which is the difference of 
ResPES spectra with the corresponding out-of-resonance spectrum. Consid- 
ering now the same intensity on Ti^+ states (Fig. 5.2l[ b)), the relative ratio 
of 5 u.c. over 3 u.c. in-gap states becomes even larger (5:1). The intensity 
of Ti^+ states marks thus the transition between the insulating and the con- 
ducting heterostructure. The resonating (red line) spectra of Fig. 5.21 give 
also a fine way to adjust the relative BE scale, by pinpointing both the 3 u.c. 
and the 5 u.c. at the maximum of Ti"^"*" states; this procedure is equivalent to 
fix the energy scale to the bulk STO Ti levels, which gives the highest spec- 
tral contribution to Ti"^"*" states. The red vertical dashed line of Fig. 5.2lK a) 
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Figure 5.19: Second ResPES data-set at L-edge on 3 u.c. LAO-STO sample. In 
panel (a) ResPES data in a 2D picture and in panel (b) a detailed view on 

Ti3+ 

defects; on the left axis XAS (black line) is depicted. In panel (c), spectra taken 
at out-of resonance (blue line, blue dashed horizontal line of panel (a)), maximum 
of Ti^"*" spectra (red line) and maximum of Ti^"*" (black line) conditions. In panel 
(d) the detail of in-gap state spectra. 



marks the reference used for the alignment, which results in the same BE 
also for the STO shallow core levels (O 2s and Sr 4p, vertical black dashed 
line). With this alignment procedure, the in-gap states of both samples are 
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Figure 5.20: CIS spectra of the second ResPES data-set; green and black solid 
lines are XAS taken at Ti L-edge on the 5 u.c. sample and on a LaTiOs crystal 
(adapted from Ref. *'^); cross markers are for the 5 u.c. CIS and triangles for the 
3 u.c; scaling factors have been used to visualize Ti^~^ CIS on the same scale of 
valence band resonances. 



found at the same energy, i.e. crossing the Fermi edge. 

High resolution spectra of Ti'^"^ states are shown in Fig. 5.22 along with 
a comparison with DFT+U calculations. To better prove the energy scale, 
in Fig. 5.22[ a) the measured Fermi edge of metallic gold deposited on the 
sample-holder at the same photon energy is given. The spectral shapes of 



116 



The LaAlOs-SrTiOs heterostructure 




Binding Energy (eV) Binding Energy (eV) 

Figure 5.21: Comparison of photoemission spectra taken in on-resonance (blue) 
and out-of-resonance (black) conditions on the 3 (thick line) and 5 (thin line) u.c. 
LAO-STO samples. The green spectra are the difference of the two and thus show 
the resonating contribution only. Vertical dashed lines mark the references that 
have been used for matching the relative energy scale. In panel B, a detail of Ti'^"^ 
defects, normalized to the VB. 



Ti^"*" states are rather similar and shows an asymmetrical profile that can be 
found in other 3d^ (for example, the titanium oxychloride, TiClO''^) highly 
correlated compound. However, in this case the in-gap state shape is due to 
the filling of empty d-states, as can be seen aligning our data with DFT+U 
calculations, shown in Fig. 5.22[ b). The DFT calculation has been carried 
out on a cubic single SrTiOa cell; the Hubbard U (U=8.0 eV) correction has 
been added to recover the correct 3.2 eV STO band-gap. The Ti^+ state 
is found exactly at 3.2 eV with respect to the STO VB-band maximum. 
The non-correlated nature of these electronic states can be further proved by 
comparison with the correlated titanium 3d states measured with ResPES in 
Nb doped SrTiOg'^^j 



ResPES data at lanthanum M4_5 edge are shown in Fig. 5.23, A strong 
resonant enhancement can be seen on the shallow core levels (La 5p); how- 
ever, the in-gap state is not resonating, indicating that La states are not 
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Figure 5.22: High resolution spectra on in-gap state taken at hi/ = 459.6 eV of 3 
u.c. and 5 u.c. LAO-STO, compared with LDA+U calculations. The calculation 



parameters are the same of Fig. 5.10 In panel (A), the Fermi edge (yellow circle) 



of metallic Au deposited on the sample-holder, taken as an energy reference. 



contributing to these electronic states and, ultimately to the transport prop- 
erties. Finally, more information can be extracted from La ResPES data by 
considering the difference between resonant and off-resonant spectra. As in 
the case of Ti 3d states, the resonating spectral wight at lanthanum edge can 
be predicted with DFT+U La 4f projected DOS. By aligning the experimen- 
tal data with calculation (in the 5 u.c. sample, peak A for titanium and peak 
B for lanthanum of Fig. 5.24), it is possible to find the mismatch of the top 
of valence bands (marked with C in Fig. 5.24) of STO and LAO separately, 
that is the band-bending. The relative shift between the 3 u.c. and the 5 
u.c. is again very low (0.15 eV) and it's comparable to the value extracted 
from XPS; this value is not enough to explain the difference in conductivity. 
However, the total shift between LAO and STO, calculated from data of Fig. 



5.24, can be quantified as 0.6 eV, in agreement with some theoretical model 
given in the literature'^®'. 

Finally, ResPES at Ti edge have been collected either in vertical (com- 
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Figure 5.23: In panel A, ResPES data and in panel B, XAS spectrum taken at 
La M4^5 edge on the 5 u.c. sample. 
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Figure 5.24: Comparison of ResPES at Ti L and La M edges on the 5 u.c. 
sample. Top four spectra (black and red lines) are experimental data take at on 
and off-resonance condition; middle (blue and green lines) spectra are the total 
DOS calculated with DFT; lower spectra (yellow and grey) show the difference 
of ResPES spectra as compared to DFT projected DOS. The DFT calculation 



parameters are the same of Fig. 5.10 
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pletely in-plane) and horizontal (partly out-of-plane) linear polarized X-rays 



see Fig. 2.6 ); the intensity ratio of Ti^"*" resonances versus the total VB spec- 



tral weight is higher with the in-plane polarization geometry (see Fig. 5.25), 



therefore the data shown on this section have been collected with that config- 
uration. However, when subtracting the off-resonance contribution from the 
on-resonance spectra, the results are virtually identical in both polarizations 



(panel B of Fig. 5.25 ). Apart from giving the optimal experimental condition 
for data collection, a polarization-dependent photoemission study can't give 
additional information concerning the physics of this system. 
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Figure 5.25: A comparison of ResPES spectra taken with different polarization. 
In panel (a), spectra taken at maximum of Ti^"'' resonance in the 5 u.c. and 3 u.c. 
LAO-STO; in the inset, a detail on the in-gap level; in panel (b), the results of 
difference of on-resonant and off-resonant spectra. 



Summarizing, the ResPES and the XAS data suggest the following: 

• There is a small but detectable Ti^"*" spectral weight at Fermi level due 
to a filling of empty STO d-bands, both in the insulating and in the 
conductive samples. Those states are measurable only in resonance 
condition; 

• The Ti'^"'" peak has a similar shape in the two sample but differs in 
intensity; 

• The CIS spectra agree with the Ti^+ picture and suggest the presence 
of intermixing; 
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• The difference in band-bending effect is virtually absent; 

• A small amount of lattice distortion has to be expected as measured 
by XLD. 

AR-XPS measurements, presented in the next paragraph, have been car- 
ried out in order to better investigate the possibility of intermixing at the 
junction. Experimental results will be discussed properly in the final section 
of this Chapter. 



5.5 AR-XPS measurements 

A cationic interdiffusion at LAO-STO junction has already been observed 
by means of X-ray diffraction, as shown in Ref. ''^ and®. In short, the 
experimental results usually display a reduction of the lanthanum content 
in the uppermost LAO layer, a strong preferential diffusion of La ions inside 



STO (Fig. 5.26 ) and a fractional occupancy of the interface Ti sites, although 
a precise quantification of the latter cannot be done. 

Given the strong angular dependence of XPS probing depth as a function 
of the emission angle, an angle-resolved XPS experiment has been carried out 
in order to reveal the chemical composition of each layer near the interface. 
The usual XPS experimental error in element quantification is usually large, 
around 10%; using a set of angular- resolved measurement it is possible to 
reduce this error, although a certain knowledge of the structure has to be 
known in advance. The 3 u.c. and 5 u.c. LAO-STO samples have been 
analyzed by angle resolved XPS, with a non-monochromatized AI-Kq X-ray 
source; the Al anode has been chosen upon Mg because of the higher photon 
energy and thus because of the higher probing depth. This kind of source 
is not focussed, thus the probing area is determined by the analyzer focus 
dimension, which is around 1 mm; since the samples had a 5x5 mm size, XPS 
signal from sample-holder hasn't been detected even at grazing emission. 
The electronic lenses of the analyzer (SCIENTA R3000 XPS/UPS/ARPES) 
allowed the selection of a small angular acceptance, that has been set to 



nearly ±1°. The experimental geometry is shown in Fig. 5.27 Both the 
analyzer and the source has been kept at fixed position, their relative angle 
a = 55° being chosen to eliminate the asymmetry correction. The alignment 
of the rotation stage has been checked by looking at the reflected spot of a 
laser beam coUimated with the analyzer nose axis. 
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Figure 5.26: Soft X-ray diffraction results on a 5 u.c. LAO-STO sample. In 
(a) the cationic occupancies and in (b) the cumulative displacement respect to 
STO layer spacing. The vertical dashed line mark the ideal interface. Taken from 
Ref.ei 




Figure 5.27: Experimental geometry for AR-XPS measurement. The analyzer- 
to-source angle a is fixed (55°) , the sample rotation is defined by the tilt angle 

e. 
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The electronic levels selected for this measurement are the Sr 3d, La 4d, 



Al 2s core levels, because of the similar binding energy (see Fig. 5.4); the 
analyzer response should be similar in such a small energy range, and the 
spectra can be collected in a single acquisition. The photoemission at Ti 3s 
peaks, which have a similar BE, has a very low cross section, so Ti data have 
been collected by considering the Ti 2p3/2 peak. When computing the total 
Ti 2p intensity, the satellite peaks have also been considered. The peak area 
has been evaluated after the removal of the satellite lines (due to the lack of 
monochromatization in the x-ray source) and after a Shirley-type background 
subtraction. 

The LAO-STO interface has been modeled with a layered structure, which 
consider each cationic plane (so the half of a layer) as a different emitter, in 
order to recover the correct La/Al and Sr/Ti ratio; however, the local value 
of IMFP is considered constant inside each oxide (LAO or STO). In the first 
approximation, the XPS intensity of a peak at a given angle is thus a sum 
of the contribution of each layer: 

I^{Ek,e) = Xm{Ek)J2ai P{\,z,9)dz (5.2) 

where m is the index of the XPS peaks, is the electron kinetic energy, 
Ui is the cationic occupancy of the half-layer i, Xm{Ek) is the cross section 



and P is the defined in Eq,2.4 The index i is running on even (right) index 
for Sr and La (Ti and Al). The half-layer thickness Zi+i — Zi is fixed to 
be equal to the lattice parameter of LAO or STO respectively, according 



to the dept z with respect to the ideal interface. The integral of Eq j5.2| 
can be performed analytically when taking A as the IMFP; the summation 
has been done on up to 50 half-unit cell. An additional correction for the 
surface environmental contamination (given by a multiplicative factor equal 
to a P{Xcont, d, 6) function) is added at every peak intensity by considering a 
fix thickness of graphitic carbon. The values of the inelastic mean free path 
have been calculated with the TPP2M® formula and are shown in the table 

o 

In a perfectly terminated LAO-STO interface the Xsto values for Al and 
La would not be needed, since photoelectrons from La and Al would travel 
only in the LAO layer; however, in the intermixing case this value should 
to be considered (e.g. for electrons emitted from La ions that have diffused 
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XPS peak 


Efc(eV) 


Alao(A) 


XSTO (A) 




La 4d 


1375 


25.09 


26.59 


36.77 


Sr 3d 


1345 


24.67 


26.14 


36.14 


Ti 2p 


1025 


20.05 


21.20 


29.25 


Al 2s 


1361 


24.90 


26.38 


36.47 



Table 5.1: List of IMFP used in AR-XPS calculations, obtained by TPP-2M 
formula. 



deeply into the STO). The contamination correction has been added to ad- 
just the angular slope of AR-XPS data; it should affect La, Sr an Al peaks in 
the same way, because of the similar kinetic energy (and thus to the similar 
Xcont)- However, since the kinetic energy of Ti 2p electrons is lower, a thick 
contamination layer can reduce the intensity ratio of Ti area versus the other 
elements. Finally, a common scaling factor is applied to match the experi- 
mental data. The experimental results, depicted with the model results in 



the case of an abrupt interface, are given in Fig. 5.28 



The peaks area is shown "as it is", without the cross-section normaliza- 
tion. The AR-XPS data clearly show the signature of X-ray photoelectron 
diffraction (XPD), consistent with the data reported in literature''^; in par- 
ticular, the XPD peaks at 0° and 45° can be detected, due to the forward- 
scattering events that occur in a cubic lattice with this kinetic energy. This 
signal modulation can increase considerably the amount of error in a fitting 
procedure. In any case, the result shown in Fig. 5.28 clearly doesn't match 
the experimental results. The calculations systematically overestimate the 
Ti and Sr contribution. Besides, an exchange between La and Sr ions at the 
interface would even increase the discrepancy with measured data (see the 
arrows in Fig. 5.28). In order to improve our model, two things have to be 
done: 



• On the experimental data side, it is needed to remove the photoelectron 
diffraction contribution from the angular dependence of the peak area; 

• On the theoretical side, it is needed to work with a better approxima- 
tions for A than the simple IMFP. 

The first issue can be solved by fitting the experimental data with a 
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Figure 5.28: Experimental AR-XPS measurement (filled circle) on the 3 u.c. and 
5.U. LAO-STO, with the calculation results (thick line). The calculation has been 
done on a perfectly abrupt LAO-STO configuration, with the implementation of 
IMFP only. The model has been normalized to La 4d data. The arrows show the 
effect of a possible La-Sr intermixing. 



function that interpolates the angular dependence, averaging out the FED 



effects; following Eq,2.4, this function can be expressed in this form: 



f{9) = Aexp- 



(5.3) 



with A and B being the parameters to be determined. The result of the data 
interpolation could in principle be fitted directly with the AR-XPS model; 
however, the estimated error on the A parameter is rather large as a results 
of the XPD modulation on the data. To better illustrate this concept, the 
experimental data on the next figures are replaced by the interpolation results 



with the Eq. 5.3 with an error bar given by: 

ferr{0) = ia^exp-^ (5.4) 



B 



where a a is the estimated error of the A parameter in the fit. 

As explained in Chapter 2, the theoretical model can be improved by 
replacing the IMFP with effective attenuation length (EAL), evaluated for 
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Figure 5.29: Calculated Effective Attenuation Length (EAL) for the measured 
XPS peaks in LAO and STO. From the NIST EAL database software. 



a thin overlayer geometry. Calculated EAL's for the LAO-STO case are 



shown in Fig. 5.29 Since Xeal is a function of the depth z, the integration 



of Eq. |5.2| has to be done numerically. As shown in Fig. |5.29[ in general 
EAL increases strongly at the surface, reflecting the fact that electrons from 
bulk STO should undergo to more elastic scattering events. The surface 
contamination has also to be evaluated trough a different approach than the 
simple graphitic layer; the carbon contamination can be estimated with an 
average EAL, that can be approximated with the phenomenological formula: 



Xcx>nt{Ek 



0.0129E^-^^^^(nm) 



(5.5) 



The parameters of Eq. 5^ have been obtained through a fitting procedure 
over many measurements on different samples 

The results of the improved model and of the data fitting are shown 



in Fig. |5.30[ A 1.2 nm coverage of hydrocarbon contamination has been 
estimated by fitting AR-XPS data, for both the 3 u.c. and the 5 u.c. LAO- 
STO samples. A much better agreement with experimental results is found, 
even without intermixing effect. In particular, the Al 2s and La 4d areas 
are in good agreement with the theoretical model, while Sr and Ti peaks 
are sometimes a little off the experimental error bar. These results could 
be considered as a proof of the effectiveness of the EAL approach. It isn't 
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possible to make a reliable structural predictions starting from the results of 
Fig. 5.30 In fact, it should be noted that there's a common multiplicative 
factor; as an example of the possible mistakes, an excess of Sr can be fitted 
either by increasing the Sr occupancy or by decreasing the occupancies of 
all the other elements. This, together with the XPD fluctuation, does not 
allow a precise fitting procedure in this case. In particular, in this case AR- 
XPS can't give a precise prediction of the depth of the intermixing, but can 
anyway discriminate whether an intermixing is present. 

These AR-XPS data have been used to confirm or refine the results ob- 
tained with other techniques. In particular, the SXRD measurements show 
a reduction of the topmost LaAlOs layer into a reduced La(i_2,)A103 and dis- 
play a cationic intermixing occurring in the firs layer below the interface, with 
the formation of a La doped STO compound. The composition proposed for 
the mixed La-Sr layer in Ref. can be expressed as Lax^Ti_3/2x] this for- 
mulation preserves the charge of a single plane, although it doesn't describe 
a complete filling of each layer. This kind of intermixing thus corresponds to 
a lattice rich of cationic vacancies. 



A refined solution for the 5 u.c. and 3 u.c. is given in Fig. 5.31 and in 
Fig. 5.32 In both samples, the top layer has been set to a partial occupancy 
x=0.8 (5 u.c.) and x=0.5 (3 u.c), the intermixing layer has been set to a 
Lao.sSro.a. With these values, the total amount of La as compared to Al is 
nearly conserved, with a reduced content of Sr on the STO termination. In 
order to adjust the Ti 2p intensity, in the 5 u.c. the occupancy axi at the 
interface has to be reduced; the best fit is with aTi=0.3, but every value 
below 0.6 is within the experimental error. Remarkably, in the 3 u.c. case, a 
very good matching is found even without reducing the Ti occupancy at the 
interface. The different Ti stoichiometry in conducting and insulating LAO- 
STO can thus be reflected on the intensity of Ti^+ states as measured by 
ResPES and XPS; however, further investigation would be needed to reduce 
the experimental uncertainty on Ti occupation. 



Finally, in Fig. |5.33| AR-XPS data taken on the insulating 5 u.c. LAO- 
STO (grown with a 10~^ mBar O2 partial pressure), compared to the con- 
ductive 5 u.c, are shown. With this photon energy (Al k^ X-ray source), the 
photoelectron diffraction peaks can be interpreted with forward-scattering 
events and thus marked with the corresponding lattice direction, as shown in 
Fig. 5.33 The possibility to observe the XPD modulation can be regarded 
as a proof of the crystalline order of these heterostructures. Consistently 
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Figure 5.30: Experimental AR-XPS measurement (filled circle) on the 3 u.c. 



and 5.U. LAO-STO, fitted with Eq. 5.3, with the calculation results (thick line). 
The confidence bands on the experimental data fitting are also shown. For the 
sake of clarity, the spectra for the perovskite A and B cations are displayed in 
different graphs. These calculations have been done on a perfectly abrupt LAO- 
STO configurations, with the introduction of the EAL. The layer occupation for 
each case is shown in the lower panel. 



5.5 AR-XPS measurements 



129 



w 
a. 
>< 

DC 
< 



•t"' ' ' ' ' ' ' ' ' I ' ' ' ' ' ' ' ' ' I ' ' ' ' ' I I I 



-1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



■ Ti 2p 

■ AI2s 




20 40 60 

Tilt angle (°) 



Tilt angle (°) 



Ti 



Sr 



in 0.8 

>, 0.6 

" 0.4 - 

-' 0.2 >— 



n 



A! 



La 

-■-1- 



10 20 30 

Half-layer number (10 is the interface) 



40 



Figure 5.31: Experimental AR-XPS measurement (filled circle) on the 5 u.c. 
LAO-STO, with the calculation results (thick line). The confidence bands on the 
experimental data are also shown. The cationic occupancy used in the calculation 
is shown in the lower panel. 
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Figure 5.32: Experimental AR-XPS measurement (filled circle) on the 3 u.c. 
LAO-STO, with the calculation results (thick line). The confidence bands on the 
experimental data are also shown. The cationic occupancy used in the calculation 
is shown in the lower panel. Note that the Ti occupancy at the interface layer is 
fuU. 
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Figure 5.33: AR-XPS data taken on conductive and insulating 5 u.c. samples. 
The lattice directions associated to the XPD forward scattering process are also 
shown. Data of the two samples are normalized to match the La 4d AR-XPS 
spectra. 



with results shown in Fig. 5.4, the Al 2s peak intensity is weaker in the in- 
sulating sample than in the conductive one, because of the different growth 
regime. By comparing the intensities of each cationic species, it seems that 
the insulating 5 u.c. sample is thus aluminium poor rather than lanthanum 
rich. The LAO-STO samples grown in 10~^ mBar (or higher) O2 pressure 
thus can't be compared with the ones grown at lower partial pressure, which 
also displays the physical properties of interest. 

Summarizing the AR-XPS results, a clear signature of La intermixing in 
the bulk STO is found in both the 3 u.c. and 5 u.c. samples, in a similar way. 
A good data matching can be found following the structural model extracted 
from SXRD measurements, which should be accompanied by a reduction of 
Ti at the interface in the 5 u.c. case. In any case, since an abrupt model does 
not match exactly with the experimental data, every physical description of 
the conductivity in the LAO-STO heterostructures should account for an 
amount of disorder at the interface. 
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5.6 Conclusions 

The most interesting signatures of the conductivity shown in this Chapter 
are the in-gap states detected with ResPES at the Ti L2,3 edge. We provided 
the first experimental evidence''^ of these states, probed with soft X-rays 
spectroscopies. FoUowing our investigation, a similar ResPES analysis has 
been carried out by A. Koitzsch et al. ® on LAO-STO interfaces grown 
with different oxygen partial pressure. They reported the presence of in-gap 
states, especially in their "oxygen poor" sample. Also in this case the oxygen 
stoichiometry seems to play a crucial role to determine the presence of Ti 
3d^ electronic levels . 

In our samples, the different spectral weight of in-gap states seem to be 
directly related to the interface conductivity. Surprisingly, the Ti^"*" states 
have also been clearly observed in the insulating sample, with the same peak 
shape and at the same energy of the conductive one. In a common metal, the 
rise of conductivity should be marked by an increasing DOS at Fermi level, 
as measured in LAO-STO, induced by a shift of the Fermi- level, which has 
not be observed in this case. The XPS binding energies are almost the same 
in insulating and conductive LAO-STO. The presence of Vq due to X-ray 
damages should be excluded, since the X-ray penetration depth is much larger 
than the LAO capping and thus a similar amount of X-ray defects should has 
been detected in each sample. Moreover, the photoemission spectra did not 
change upon the exposition to the synchrotron X-ray beam over the entire 
period of the data acquisition (nearly 5 days). 

As a starting point we should consider that both samples are correlated 
systems, with a much higher effective mass in the "insulating" sample; in that 
case, however, a different quasiparticle peak-shape at Fermi level should be 
observed, as can be seen in 3d^ correlated metals (for example, in V2O3®) 
and in theoretical models (for example, the Dynamical Mean Field The- 
ory®!'). Moreover, as shown with DFT+U DOS calculations, the in-gap 
states can be described with a simple band-filling of the empty STO lev- 
els®. 

Since the difference in the Ti^"*" states of 3 and 5 u.c. samples is only 
a scaling factor without a detectable energy shift, the conductivity may be 
related to a phase separation mechanism which, in turn, is driven by the 
thickness of the LAO layer. In any case, the problem is to find a suitable 
"doping" mechanism that can justify the presence of these conductive Ti^"*" 
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Figure 5.34: Doping versus temperature phase diagram of SrTiOs/LaAlOs in- 
terfaces. Tlie doping scale consists of tliree possible contributions: intrinsic carrier 
doping from the electronic reconstruction (xj„t), carrier doping from oxygen vacan- 
cies (xoa ) and carrier doping by applied electric fields (xjjeid)- Observed transport 
effects are insulation at p-type interfaces ''^(xjO), 2D interface metallicity , 3D 
bulk metallicity*tZ3j negative magneto resistance''^ below 70 K, magnetic hys- 
teresis''^ below 300 mK, and superconductivity''^ below Tc=200 mK. The ex- 
act position of Xiut in the phase diagram is not yet known and can vary due to 
other polarity discontinuity compensating mechanisms (such as xoj ). Taken from 
Ref.CD. 



states at the LAO-STO interfaces. 

The different conductivity mechanism found in the hterature have been 



summarized''^ by M. Hujiben et ah in the phase diagram of Fig. 5.34 Three 
different type of doping are underhned: the intrinsic doping (xj„t), related to 
the electronic reconstruction, the extrinsic doping (xq^, mostly due to Vo) 
and the external electric field {x field). It should be noted that an electrical- 
field driven transition from the superconductive to the metallic state has 
been reported thus justifying the presence of the xjjgZd parameter in Fig. 



5.34[ The 3D metallic phase correspond to a high density of Vo, due to 
the low P02 during the growth. The superconductivity has been observed 
for an intermediate P02 regime , while in the optimal regime the localized 
magnetic moments have been detected. The sample analyzed in this thesis 



work belong to the latter group. The value of Xi^t in Fig. 5.34 is still 
unknown; in particular, it cannot be put in the intermediate region because 
oxygen vacancies can have a non negligible role in the superconductivity. 

The observed in-gap states could be ascribed for instance to a cationic dis- 
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order at the interface; in fact, STO can became conductive even for small La 
doping. In order to justify the 4 u.c. thickness threshold, the disorder should 
be dependent to the number of LAO capping layer; however, the cationic 
interdiffusion seems to be similar in conductive and insulating samples, as 
shown in section 5.5 and in the literature by N. Nakagawa et al. ®. 

The metallic areas can also be generated by the presence of residual oxy- 
gen vacancies near the interface; as in the case of Ti02-5, an oxygen vacancy 
could create localized defect states, that above certain concentrations can 
merge to yield a metallic band. In a recent article''^, Z. Zhong and co- 
workers report DFT calculations suggesting that the formation energy of 
oxygen vacancies in LAO-STO interfaces is thickness dependent. However, 
our samples should also display the presence of isolated magnetic moments 
that in turn cannot be provided by isolated oxygen vacancies alone, similarly 
to the Ti02_5 case of Chapter 4. The magnetism could appear only in the 
case of clustered defects at the interface''^. The possibility to reduce the 
influence of Vo by mean of the growth conditions and the sharp thickness 
threshold suggest that interface Vo can play a minor role in the conductivity, 
at least in our samples. 

Finally, in the electronic reconstruction scenario the relaxation of the 
polar discontinuity is achieved by a 1/2 e~ charge-transfer from the surface 
to the interface, which in fact should result in the presence of reduced Ti'^'^+ 
in the uppermost STO layer "l^. A practical implementation of this concept 
is the band-bending effect ® , which has never been observed in these systems 
by means of X-ray spectroscopies''^. Moreover, ab-initio calculations carried 
out on a LAO-STO interface with different types of cationic disorder''^ show 
that the intermixing can be thermodynamically favorable and can reduce the 
LAO polarity effects. 

Another kind of electronic reconstruction is a ferroelectric-like distortion 
of the first LAO layer at the interface. As expected, the average c cell pa- 
rameter in the capping LAO (3.79 A) is usually lower than the STO planar 
distance''^ (3.89 A); however, SXRD and HRTEM'^S measurements show 
that the first unit cell at the interface is elongated by 4%-9% (up to 4.1 A). 
This elongation is explained by the authors with the presence of Ti'^+ ions 
at the interface (thus with a larger ionic radius with respect to Ti^^), in 
a sort of Jahn- Teller effect. A structural deformation due to the strain in 
LAO capping that is spread also in STO has also been reported in a recent 
investigation with high-resolution XRD®. As already pointed out in Sec- 
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tion 5.3, the distortion of tlie oxygen octaliedral cage around tlie interface 
Ti atoms has already been measured with XLD'ED and in a similar way in 
our samples. Surprisingly, Ti^+ electronic levels cannot be observed by XAS, 
and XLD results can be explained completely with Ti^"*" CT-multiplet calcu- 
lation only. However, CIS spectra clearly agrees with a Ti^"*" picture of the 
conductive Ti 3d states. In any case, in the electronic reconstruction frame 
the presence of the Ti^"*" states at Fermi edge in insulating LAO-STO cannot 
be explained, at least without assuming an inhomogeneous distribution of 
the distorted sublattice at the interface. 

Another interesting theory that could fit our experimental results has 
been proposed by N. C. Bristowe and co-workers®; by recognizing that 
most of the measurement on this oxides have been done at atmospheric pres- 
sure (or, as in this case, with an environmental atmosphere surface contam- 
ination), the doping at the LAO-STO interface could be given by a surface 
contamination or surface defects. 

The relaxation of the polar discontinuity can thus be achieved by the same 
charge-transfer process of the electronic reconstruction model, but without 
the need of a band bending, which should also imply a conductive surface. In 
this model, the oxygen reduction at the surface should create a potential well 
for the induced localized charge in the 3d-level at the interface; in the authors 
view, the localization radius would be inversely proportional to the depth of 
the well, i.e. to the LAO thickness: hence the reason of the conductivity. A 



pictorial representation is given in Fig 5.35 



This model can also explain a peculiar effect that has been reported on 
insulating LAO-STO interface. In fact, Y. Xie et al. have shown that it 
is possible to draw (and erase) conductive lines on LAO-STO by moving a 
biased AFM tip on the LAO surface. The authors suggest that the driving 
mechanism of this process is the creation of defects at the surface; in the 
surface-reduction model, the biased tip just raises (or decreases, according 
to the bias sign) the vacancy stability. Recently published ab-initio DFT cal- 
culations of the Vo formation energies in the LAO-STO further support this 
description®. A possible way to check the influence of the surface reduction 
on the LAO-STO conductivity could be to use LAO-STO junction as a "gas 
sensor" ; it should be possible to measure a variable sample resistance upon 
the exposition to different gases or to directly measure the Ti-3d spectral 
weight with the aid of ResPES in a synchrotron facility. 

In conclusion, the LAO-STO system displays a wealth of physical and 
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The LaAlOs-SrTiOs heterostructure 
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Figure 5.35: Schematic view of the thickness dependence on interface electrons 
locahzation from the model of Ref . . Trapping potential V created by a surface 
O vacancy as seen by interface electrons versus distance within the interface plane 
X (x=0 is directly below the vacancy) for film thickness d=a (deepest), 2a, 3a, 
4a (shallowest), with harmonic estimates of corresponding donor ground states 
[taking me// = 3me]. Inset: Sketch of range and density of trapped states. 

chemical properties, which challenge both theoretical and experimental physi- 
cist. A single interface "doping" source is not enough to explain the various 
phenomena that occur at this interface. In fact, the theoretical models that 
involve perfectly abrupt interfaces and clean surfaces cannot completely ex- 
plain the different experimental evidences. In this scenario, surface effects 
and electronic phase separation are thought to play a major role, so far often 
disregarded. 
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